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I.  INTRODUCTION 


A  review  of  recent  observational  studies  of  the  main 
phenomena  of  this  study,  the  Mei-Yu,  and  the  associated 
disturbances  will  be  presented  first.  There  follows  a 
review  of  the  pertinent  theoretical  studies.  A  description 
of  the  objectives  of  this  numerical  study  of  the  Mei-Yu 
front  will  then  be  given. 

A.  OBSERVATIONAL  STUDIES 
1.  Mei-Yu 

Over  eastern  Asia,  the  transition  period  from  the 
northeast  winter  monsoon  to  the  southwest  summer  monsoon 
possesses  a  distinctive  climatological  phenomenon  which  is 
important  for  regional  weather  forecasting.  The  seasonal 
rainfall  reaches  a  maximum  during  this  period  (Fig.  1.1). 
The  rainfall  during  this  period  is  called  "Mei-Yu"  (plum 
rain)  in  China  and  "Baiu"  in  Japan.  During  this  period,  the 
rainfall  may  be  continuous  or  intermittent  for  several  days 
to  a  few  weeks  and  is  characterized  by  frequent  rain-showers 
and  thunderstorms,  with  rainfall  rates  of  up  to  a  few 
hunc'-ed  millimeters  per  day  (e.g.,  Chen,  1983;  Ninomiya, 
1980).  The  Mei-Yu  season  first  appears  in  southern  China  in 
May,  and  influences  the  weather  over  southern  China  and 
Taiwan  from  mid-May  to  mid-June  (Chen  and  Chi,  1980).  It 
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then  migrates  over  the  Yangtze  Valley  and  Japan  from  mid- 
June  to  mid-July. 

Synoptical ly ,  the  Mei-Yu  rainfall  is  associated  with 
the  repeated  occurrence  of  a  front  that  develops  in  the 
midlatitudes  and  slowly  moves  southeastward  to  a  quasi¬ 
stationary  position  extending  from  southern  Japan  to 
southern  China.  Satellite  .pictures  reveal  a  nearly 

continuous  cloud  band  accompanying  this  front,  which  is 
often  called  a  Mei-Yu  front  or  a  Baiu  front.  The  frontal 
position  coincides  with  the  early  summer  monsoon  trough  in 
eastern  Asia  and  its  existence  may  be  viewed  as  a 
manifestation  of  the  trough.  The  beginning  of  the  rainy 
season  also  signifies  the  onset  of  the  northern  summer 
monsoon  in  this  part  of  the  world.  The  Mei-Yu  often 
undergoes  sequential  jumps  between  May  and  July  (Fig.  1.2; 
Tao  and  Chen,  1987)  depending  upon  the  strength  and  the 
position  of  the  North  Pacific  Ocean  subtropical  high  (Tao 
and  Ding,  1981).  As  the  Mei-Yu  onset  occurs  in  tiie  Yangtze 
River  region  in  June,  the  Indian  monsoon  is  also  being 
established  and  the  westerly  jet  stream  begins  to  retreat  to 
the  north  of  the  Tibetan  Plateau.  At  the  same  time,  the 
subtropical  high  in  the  western  Pacific  area  moves  northward 
to  2  0  -  2  5  °  N  .  The  Mei-Yu  season  lasts  through  inid-July,  when 

the  subtropical  high  moves  northward  to  30°N.  At  that  time, 
the  rain  belt  is  to  the  north  of  the  Yellow  River,  which 
signifies  the  beginning  of  heavy  rainfall  period  in  northern 


China  and  the  ending  of  the  Mei-Yu  in  subtropical  China  (Tao 
and  Chen ,  1987). 

2 .  Mei-Yu  Front 

The  climatology  of  the  Mei-Yu  front  over  southern 
China  and  the  adjacent  areas  was  examined  by  Chen  and  Chi 
(1980)  for  the  periods  of  15  May-15  June,  1968-1977.  They 
found  that  the  Mei-Yu  f rontogenes  is  takes  place  in  the 
subtropical  part  of  southern  China  to  the  south  of  ~35°N. 
All  the  Mei-Yu  fronts  were  formed  in  the  area  of  20°-35°N 
and  100°-130°£  (Fig.  1.3).  The  lifetime  of  individual  cases 


of  Mei-Yu  fronts  ranges 

from  3  days 

to  22 

days 

with 

an 

average  of  8 

days  . 

The  kinematic 

background 

for 

the 

f  rontogenes  is 

is  the 

deformation 

field 

between 

the 

subtropical  high  over  the  Pacific  and  the  migratory  high 
over  central  China.  They  suggested  that  the  sea-surface 
temperature  contrast  in  the  nearby  South  and  East  China  Seas 
may  also  be  a  favorable  factor  for  the  f ron t ogenes i s . 

In  contrast  to  the  polar  front,  the  temperature 
contrast  across  the  Mei-Yu  front  is  typically  small,  but 
there  is  always  a  marked  moisture  gradient  near  the  frontal 
zone  (Akiyama,  1973a;  Chen  and  Tsay,  1978;  Matsumoto  et  al., 
1970,  1971;  Tao,  1980).  An  analysis  of  the  kinematic 
structure  of  the  Mei-Yu  frontal  system  over  Southeast  Asia 
was  carried  out  by  Chen  (1978).  The  Mei-Yu  frontal  system 
possesses  a  remarkable  shear  line  which  roughly  coincides 
with  a  trough  at  850  and  700  mb,  although  the  associated 
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geopotential  gradient  is  rather  weak  in  the  immediate 
vicinity  of  trough.  The  shear  line  approximately  coincides 
with  the  maximum  north-south  gradients  of  mixing  ratio  and 
equivalent  potential  temperature.  To  the  south  of  the  shear 
line,  a  low  level  jet  (LLJ)  that  is  characterized  by  maximum 
mixing  ratio  and  high  equivalent  potential  temperature  is 
found  near  700  mb  (Fig.  1.4).  .  Strong  cyclonic  vorticity, 
horizontal  convergence  and  upward  motion  are  observed  along 
and  to  the  south  of  the  shear  line.  The  cloud  band 
approximately  coincides  with  the  area  of  maximum  cyclonic 
vorticity,  horizontal  convergence  and  upward  motion. 
Intense  convective  clouds  in  the  southern  portion  of  the 
cloud  band  tend  to  be  associated  with  the  maximum  values  of 
these  kinematic  parameters  as  well  with  relative  humidities 
greater  than  80%  in  the  lower-middle  troposphere.  Clear 
skies  over  the  East  China  Sea  and  its  vicinity  are 
accompanied  by  downward  motion  of  dry  air  (relative  humidity 
less  than  30%  at  850  mb,  less  than  10%  at  700  and  500  mb). 

Chen  and  Chang  (1980)  studied  the  structure  and 
calculated  a  vorticity  budget  for  a  Mei-Yu  system  during  the 
period  10-15  June  1975.  During  the  mature  and  decaying 
stages  of  the  monsoon  trough,  the  structure  of  the  eastern 
section  (southern  Japan)  and  the  central  section  (East  China 
Sea)  is  different  from  that  of  the  western  section 
(southeastern  China).  In  the  eastern  and  central  section, 
the  structure  resembles  a  typical  midlatitude  baroclinic 


front  with  a  strong  vertical  tilt  toward  an  upper-level  cold 
core  low  and  a  strong  horizontal  temperature  gradient.  On 
the  other  hand,  the  wester  n  section  resembles  a  tropical 
disturbance  with  an  equivalent  barotropic,  warm  core 
structure,  a  relatively  weak  horizontal  temperature  gradient 
and  a  rather  strong  horizontal  wind  shear  in  the  lower 
troposphere.  In  the  mature  stage,  the  upward  motion  reaches 
maximum  values  at  500  mb.  It  is  located  over  the  850  mb 


monsoon 

trough  in  the  western 

(  W) 

section, 

4°  latitude  south 

of 

the 

trough  in  the  central 

(  C  ) 

and  3° 
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the 

trough  in  the  eastern 

(  E  ) 
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trough,  downward  motion 

pr eva l Is 

with  a  maximum  at 

400 

mb 

in  sections  C  and  E, 

and  at  700 

mb  in  section  W. 

This  well-organized  vertical  velocity  couplet  is  indicative 
of  a  well-developed  secondary  circulation  normal  to  the 
trough.  In  the  decaying  stage,  the  circulation  becomes 
substantially  less  vigorous  in  sections  W  and  C,  but  it 
maintains  intensity  in  section  I'.  1'sav  and  Chen  <1980). 
analyzed  the  processes  which  force  the  vertical  motion  in 
sections  W,  C  and  E  .  They  found  the  eastern  portion 
partially  maintains  the  characteristics  of  a  midlatitude 
system,  while  the  western  portion  behaves  similar  to  a 
t  rop i ca 1  system. 

The  Mei-Yu  front  can  be  several  thousand  kilometers 
long,  and  therefore  it  is  a  large-scale  monsoonal 


phenomenon.  However,  the  cross-front  scale  is  considerably 


smaller,  which  is  within  the  mesoscale  range.  Among  the 
many  mesoscale  features  that  are  observed  near  the  Mei-Yu 
frontal  zone,  the  most  important  are  the  organized  mesoscale 
convective  systems  embedded  within  the  broad  Mei-Yu  cloud 
band,  which  lies  on  the  warm  side  of  the  front  (Fig.  1.5). 
These  convective  systems  tend  to  move  along  the  front  from 
west,  to  east  and  produce  local.ly  heavy  rainfall  rates  of 
>  100  mm  day-1  (Akiyama,  1974,  1978;  Chen,  1977;  Ninomiya 
and  Akiyama,  1971,  1972,  1973,  1974). 

3 .  Low  Level  Jet  (LLJ) 

To  the  south  of  the  Mei-Yu  (or  Baiu)  front,  there  is 
often  a  strong  ageostrophic  low  level  jet  (LLJ)  which  is 
usually  observed  between  1.5  and  3  km  (850-700  mb)  and  is 
strongly  correlated  to  the  heavy  rainfall  (Akiyama,  1973b; 
Chen,  1977,  1986;  Matsuraoto,  1972;  Matsumoto,  et  al.,  1971; 
Ninomiya  and  Akiyama,  1974;  Tao  and  Chen,  1987).  Akiyama 
(1973a)  observed  that  the  heavy  rainfalls  in  the  Japan  Baiu 
season  occur  most  frequently  to  the  north  of  the  low  level 
jet  stream  and  seldom  occur  in  areas  more  than  200  km  south 
of  the  jet  axis.  Matsumoto  (1973)  found  that  the  peak  of 
heavy  rainfall  precedes  by  6  h  the  peak  of  the  low  level  jet 
speed.  He  also  found  that  the  peak  of  heavy  rainfall  is  in 
phase  with  the  maximum  acceleration  of  the  LLJ. 

In  the  Great  Plains  and  Midwest  areas  of  the  United 


States,  the  low  level  jets  are  a  climatological  feature  with 
the  maximum  frequency  during  the  spring  and  summer  months 


(Bonner,  1968).  The  low  level  jets  in  this  area  are 
basically  a  boundary  layer  phenomenon  that  are  characterized 
by  a  diurnal  oscillation.  The  jet  reaches  maximum  intensity 
by  early  morning,  and  is  associated  with  a  nocturnal 
temperature  inversion  (e.g.,  Blackadar,  1957;  Bonner,  1968; 
Holton,  1967;  Lettau,  1967;  Wexler,  1961).  Different 
theories  (discussed  later  in  section  B)  have  been  proposed 
to  explain  the  generation  of  the  low  level  jet  and  the 
seasonal,  temporal  and  geographic  variations.  In  the  British 
Isles,  a  low  level  jet  in  conjunction  with  midlatitude 
cyclones  was  studied  by  Browning  and  Harrold  (1970)  using 
Doppler  radar  observations.  Browning  and  Pardoe  (1973) 
found  that  for  six  cases  of  ana-cold  fronts  occurring  in  the 
colder  months  in  British  Isles,  the  low  level  jet  had  a 
maximum  velocity  in  the  layer  900  to  850  mb  of  25  to  30  m  s1 
and  occurred  during  night  or  day. 

4 .  Mesoscale  convective  systems  (MCS’s) 

Rainfall  characteristics  over  Taiwan  and  vicinity 
were  studied  by  Chen  (  1977)  and  Chen  and  Tsay  (  1978)  using 
dense  hourly  rainfall  reports,  radar  echoes  and  satellite 
pictures.  The  low  level  southwesterl ies ,  especially  the 
LLJ,  tends  to  create  potential  instability  to  the  warm  side 
of  Mei-Yu  front.  Continued  large-scale  ascent  then  leads  to 
the  release  of  potential  instability  through  organized 
mesoscale  convective  systems.  The  mesoscale  convective 


systems  embedded  within  the  broad  Mei-Yu  cloud  band  can  be 


of  several  different'  sizes  (raeso  -a,  -  6,  or  — y).  A  case 
study  by  Chen  and  Tsay  (  1978)  revealed  the  existence  of 
these  mesosystems  in  Taiwan  with  wavelengths  of  60,  85  and 
170  km  and  the  corresponding  periods  of  2.6,  3.5  and  7.1  h. 
A  spectrum  analysis  (Tao,  1980)  of  radar  echoes  in  the 
rainbands  along  the  Mei-Yu  front  over  the  Chinese  mainland 
also  showed  the  existence  of  organized  mesoscale  convective 
systems  with  wavelengths  of  150,  84  and  50  km.  Akiyama 
(1978)  found  that  the  echoes  within  the  rainband  in  the  Baiu 
season  in  Japan  are  organized  into  mesoscale  clusters  with 
spatial  intervals  on  the  order  of  100-150  km. 

S imul at  ions  of  these  mesoscale  convective  systems 
require  a  three-dimensional  model.  In  our  two-dimensional 
simulations,  we  will  not  study  these  phenomena. 

B.  THEORETICAL  STUDIES 
1 .  Front 

Quas i-geos t rophi c  frontal  theories  (Stone,  1966; 
Williams,  1968;  Williams  and  Plotkin,  1968)  showed  that 
large-scale  deformation  fields  can  lead  to  surface 
f rontogenes is  and  result  in  strong  horizontal  temperature 
gradients  that  are  characteristic  of  fronts.  However,  the 
structure  of  the  front  and  the  associated  secondary 
circulation  do  not  tilt  with  height.  It  also  lacks  cyclonic 
wind  shear  at  the  front  and  has  several  other  unrealistic 


features,  including  the  infinite  time  required  for  the  front 


to  form.  Later  studies  using  primitive  equations  (Williams, 
1972)  and  the  semi-geostrophic  equations  (Hoskins  and 
Bretherton,  1972)  removed  these  unrealistic  features  by 
including  the  ageostrophic  advections  of  temperature  and 
momentum.  In  particular,  Hoskins  and  Bretherton  (1972)  used 
the  set  of  semi-geostrophic  equations  to  show  that  the 
nonlinear  effects  of  the  ageostrophic  advections  can  be 
represented  by  the  quas i -geos t r oph i c  solutions,  if  these 
solutions  are  written  in  the  geostrophic  coordinates  and 
then  transformed  to  the  physical  space.  Later,  the  semi- 
geostrophic  theory  was  expanded  to  three  dimensions  by 
Hoskins  (1975)  and  Hoskins  and  Draghici  (1977). 

Frictional  and  diabatic  effects  also  have  been 
studied  by  a  number  of  investigators.  Williams  (1974)  found 
that  vertical  mixing  can  determine  the  slope  of  a  front, 
and,  together  with  horizontal  mixing,  can  determine  the 
cross-width  of  the  frontal  zone.  Keyser  and  Anthes  (1982) 
showed  that  the  proper  simulation  of  the  concentrated  upward 
motion  in  the  surface  frontal  zone  depends  on  the  planetary 
boundary  layer  parameterization.  Williams  et  al.  (1981), 
Ross  and  Orlanski  (  1978),  and  Hsie  and  Anthes  (1984)  found 
that  latent  heat  release  intensifies  both  the  upper  level 
westerlies  and  the  low  level  easterlies  poleward  of  the 
surface  front,  and  the  ageostrophic  secondary  circulation 


across  the  front. 


2 .  Symmetric  Instability 

Within  extrat ropical  cyclones,  banded  structures  in 
clouds  and  precipitation  have  been  observed  both  in 
conditionally  stable  and  unstable  atmospheres  (Hobbs,  1978). 
Bennetts  and  Hoskins  (1979)  suggested  that  these  bands  are  a 
result  of  conditional  symmetric  instability.  Emanuel  (1979) 
demonstrated  that  circulations  resulting  from  this 
instability  are  fundamentally  mesoscale  in  character. 
Recent  observations  indicate  that  conditional  symmetric 
instability  is  responsible  for  some  precipitation  bands 
within  extratropical  cyclones  (Bennetts  and  Sharp,  1982; 
Moore,  1985;  Parsons  and  Hobbs,  1983;  Sanders  and  Bosart, 
1985)  . 

Bennetts  and  Hoskins  (1979)  and  Emanuel  (1982) 
assessed  the  degree  of  conditional  symmetric  instability 
based  on  infinitesimal  displacements  in  saturated 
atmospheres.  Emanuel  (1983a, b)  considered  parcel 
instability  which  leads  to  a  measure  of  the  stability  of  the 
moist  baroclir.  ic  atmosphere  to  finite  slantwise  reversible 
displacements  of  a  two-dimensional  air  parcel  or  air  "tube". 
Dynamically,  inertial,  symmetric  and  convective 
instabilities  are  very  closely  related .  Each  of  these 
instabilities  may  be  thought  of  as  resulting  from  an 
imbalance  of  body  forces  acting  on  a  fluid  parcel.  In  the 
case  of  convective  instability,  the  gravitational  force  acts 
vertically.  In  the  case  of  inertial  instability,  the 


centrifugal  force  acts  horizontally.  Symmetric  instability 
simply  represents  motion  resulting'  from  a  combination  of 
these  forces,  and  thus  the  motion  is  "slantwise."  Following 
Emanuel  (1983a, b),  it  can  be  shown  that  slantwise 
instability  requires  the  slopes  of  the  pseudo  angular 
momentum  (M)  surfaces  to  be  smaller  than  those  of  the 
equivalent  potential  temperature  (  0  e)  surfaces.  This  is 
equivalent  to  stating  that  0e  decreases  along  the  constant 
M  surfaces  upward  or  that  M  decreases  in  the  cross-front 
direction  toward  warm  air  along  a  9  f-surface.  From  the 
definition  of  M,  the  latter  implies  that  the  absolute 
vorticity  is  negative  along  6  e-surfaces. 

A  case  study  by  Emanuel  (1983b)  of  a  strongly 
baroclinic  flow  with  observed  linear  bands  of  clouds  and 
precipitation  reveals  that  the  atmosphere  is  initially 
conditionally  unstable  to  upright  convection  and  tends  to 
become  neutral  with  respect  to  slantwise  displacements. 
Such  an  atmosphere  possesses  moist  adiabatic  lapse  rates 
along  M  surfaces,  even  though  it  may  be  substantially  stable 
to  vertical  displacements.  In  this  case,  conditionally 
unstable  baroclinic  flow  undergoes  a  form  of  slantwise  moist 
adiabatic  adjustment  that  renders  the  atmosphere  neutral  to 
slantwise  reversible  displacements.  The  adjustment  uses  up 
the  total  slantwise  convective  available  potential  energy 


(SCAPE)  rather  than  just  the  convective  (gravitational  part 
of)  available  potential  energy  (CAPE)  for  upright 


convection. 


After  such  an  adjustment,  the  SCAPE  of  the 


atmosphere  is  very  close  to  zero  while  the  CAPE  actually 
becomes  negative. 

Emanuel  (1985)  carried  out  a  diagnostic  study  using  a 
semi-geostrophic  model  with  small  potential  vorticity  (i.e., 
small  symmetric  stability)  in  the  region  of  upward  motion 
and  not  so  small  potential  vorticity  in  the  region  of 
downward  motion.  He  demonstrated  that  frontal  forcing 
results  in  a  strong,  concentrated,  sloping  updraft  that  is 
located  ahead  of  the  region  of  maximum  g'eostrophic 
f r on t ogene t i ca 1  forcing.  The  distance  between  the  "front" 
and  the  updraft  will,  in  the  atmosphere,  have  a  typical 


value  on  the 

order 

of  50-200  km. 

This 

was  observed  by 

Sanders  and 

Bosart 

( 1985) , 

who 

f  ound 

that 

slantwise 

convection  can 

persist  even 

after 

a  condition 

of  near- 

neutrality  has  been  achieved  under  f r on t ogenet i ca 1  forcing. 
It  has  also  been  shown  by  Thorpe  and  Emanuel  (1985)  that  in 
an  atmosphere  which  has  small  stability  to  slantwise 
convection,  f rontogenes  is  proceeds  at  a  greatly  increased 
rate  (due  to  latent  heat  release)  and  the  horizontal  sc:  ale 
of  the  ascent  in  the  warm  air  is  considerably  reduced. 
Cross-frontal  circulations  produced  by  g'eostrophic  forcing 
are  similar  in  structure  to  slantwise  convection,  despile 


the  (small)  stability  to  the  latter. 


3.  Low  Level  Jet  (LLJ) 

In  the  vicinity  of  Japan,  the  Baiu  front  is 
accompanied  by  an  ageostrophic  LLJ  near  700  mb,  which  is 
oriented  along  the  heavy  rainfall  zone.  The  ageostrophic 
LLJ  is  suggested  to  be  a  result  of  the  downward  transport  of 
horizontal  momentum  caused  by  cumulus  convection  (Akiyama, 
1973a;  Matsumoto,  1972;  Matsumoto  and  Ninomiya,  1969; 
Matsumoto  et  a l . ,  1971;  Ninomiya,  1971;  Ninomiya  and 
Akiyama,  1974).  In  China,  several  different  relationships 
between  the  low  level  jet  and  heavy  rains  have  been  observed 
in  the  Mei-Yu  season.  In  one  of  these,  a  LLJ  is  coupled  and 
parallel  to  the  upper  tropospheric  jet.  The  formation  and 
strenthening'  of  the  LLJ  is  often  associated  with  the  heavy 
rainfall.  Chen  (1982)  suggests  that  coupling  of  upper  and 
lower  level  jets  can  be  explained  by  an  unstaole  inertia- 
gravity  wave  that  arises  from  a  thermal  wind  adjustment 
process  that  corrects  an  imbalance  due  to  the  advection  in 
the  entrance  region  of  the  upper  level  jet.  This  coupling 
is  in  contrast  to  that  studied  by  Uccellini  and  Johnson 
(1979)  in  the  Great  Plains  of  the  United  States,  where  the 
LLJ  tends  to  be  perpendicular  to  and  in  the  exit  region  of 
the  upper  level  jet. 

Several  theories  have  also  been  proposed  to  explain 
the  formation  and  diurnal  oscillation  of  the  LLJ  in  the 
Great  Plains  region  of  the  United  States.  Black adar  (1957) 
suggested  that  oscillation  in  the  boundary-layer  mixing 


processes  is  responsible  for  the  diurnal  variation  of  the 
I.LJ  intensity  which  reaches  maximum  values  by  the  early 
morning'.  The  temporal  and  geographic  preference  of  the  LLJ 
have  been  further  explained  by  diurnal  radiation  fluxes  over 
sloped  terrain  (Holton,  1967;  Hsueh,  1970;  Lettau,  1967), 
and  topographical  characteristics  (Paegle  and  Rasch,  1973; 
Paegle,  1978).  Observations  also  showed  that  low  level  jets 
may  also  develop  in  response  to  synoptic  or  subsynoptic 
processes,  particularly  through  a  response  to  leeside 
cyclogenesis  that  is  common  to  the  Great  Plains  (Djuric  and 
Damiani,  1980;  Djuric  and  Ludwig,  1983:  Naistat  and  Young, 
1973;  Newton,  1956;  1967;  Reiter,  1969;  Uccellini,  1980). 
Numerical  simulations  and  a  case  study  by  Uccellini  and 
Johnson  (1979)  suggested  that  the  development  of  the  LLJ  is 
coupled  to  the  upper  tropospheric  jet  streak  by  a  mass 
adjustment  process  within  the  exit  region  of  the  streak. 

C.  OBJECTIVE  OF  THIS  STUDY 

The  objective  of  this  study  is  to  use  a  two-dimensional 
f rontogenes is  model  similar  to  that  of  Williams  et  al  . 
(1981)  to  study  the  quasi-steady  state  structure  of  the  Mei~ 
Yu  front  and  the  LLJ.  The  model  contains  a  convective- 
adjustment  type  cumulus  heating  parameterization  and  a  bulk 
aerodynamics  formula  for  the  boundary  layer.  To  study  the 
possibility  of  cumulus  momentum  transport  or  the  diurnal 


oscillation  processes,  additional  parameterizations  would  be 
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necessary.  Since  considerable  uncertainties  exist  in  both 
the  cumulus  momentum  transport  and  the  diurnal  radiation 
parameterizat ions ,  it  is  desirable  to  carry  out  a  modeling 
study  excluding  these  processes  to  see  if  the  LLJ  may  be 
developed.  One  particular  motivation  for  doing  this  is  the 
tropical  nature  of  the  Mei-Yu  front.  The  intense  convection, 
which  is  correlated  with  the  occurrence  of  strong  low  level 
vorticity,  points  to  a  different  possibility.  The  large 
latent  heat  release  provides  an  energy  source  which  may,  in 
certain  situations,  be  converted  to  the  kinetic  energy  of 
the  LLJ.  This  possibility  is  worth  investigation  especially 
in  light  of  the  recent  symmetric  instability  studies,  which 
show  that  slantwise  convection  may  occur  in  a  frontal  region 
even  if  the  symmetric  instability  condition  is  neutral  or 
slightly  stable. 

Therefore,  the  following  topics  are  pursued  in  this 

study: 

1)  Study  of  the  structure  and  dynamics  of  a  steady 

front,  which  is  investigated  by  varying  the  vertical 
stability,  moisture  content,  Coriolis  parameter  and 
boundary- 1 ayer  fluxes  in  the  model;  and 

2)  Study  of  the  relationship  between  the  occurrence  of 

the  LLJ  and  the  convective  activity  in  a  Mei-Yu 

front.  In  particular,  the  possible  role  of  a 

secondary  circulation  equatorward  of  the  front  caused 


by  convection  will  be  explored. 


In  Chapter  II,  the  basic  formulation  of  the  model  is 
presented.  In  Chapter  III,  the  results  of  the  numerical 
experiments  are  described,  but  the  main  interpretation  of 
the  results  with  respect  to  the  dynamics  of  the  LLJ  is 
discussed  in  Chapter  IV.  Chapter  V  gives  the  summary  and 


conclusions. 


Fig.  1.1 


Mean  daily  rainfall  (y-axis,  mm)  for  the  period  of  1951-1970 
at  Taichung,  Taiwan,  which  is  a  typical  station  having  a 
Mei-Yu  rainfall  maximum  (Chen,  1977). 


Fig.  1.3 


Frequency  of  frontogenesis  at  1°  longitude  by  1°  latitude 
grid  intervals  during  the  Taiwan  Mei-Yu  period  of  May  15  to 
June  15,  1968-1977  (Chen  and  Chi,  1980). 


Vertical  cross  section  of  (a)  east-west  wind  (u;  solid)  and 
north-south  wind  (v;  dash)  components  in  m  s~l;  (b)  total 
wind  speed  (solid)  in  m  s'1  and  mixing  ratio  (dash)  in 
g  kg-1  along  120°E  at  1200  GMT  11  June  1975.  J  indices 
from  left  to  right  increase  from  south  to  north  with  a  grid 
distance  of  240  km.  The  Mei-Yu  front  at  850mb  is  near  J=7 
(Chen,  1978). 


This  study  uses  a  modified  version  of  the  model 
developed  by  Williams  (1974)  and  Williams  et  a  1 .  (1981). 
The  frontogenesis  in  the  model  is  forced  by  a  nondiver gent 
basic  wind  field  that  contains  stretching  deformation  and 
includes  an  Ekman  layer  structure.  The  model  employs 
vertical  and  horizontal  diffusions  of  potential  temperature, 
moisture  and  velocity.  The  diffusion  coefficients  vary  only 
with  elevation.  Pressure  coordinates  are  used  in  the 
vertical  and  a  tropopause  is  simulated.  This  should  lead  to 
a  more  accurate  simulation  of  the  upper  regions  where  latent 
heating  will  be  important.  The  time-dependent  quantities  in 
the  model  are  functions  of  y  (cross-front)  and  p  and  the 
model  is  integrated  until  a  quasi-steady  state  is  achieved. 
A  planetary  boundary  layer  is  added  with  aerodynamic  bulk 
formulae  for  the  surface  fluxes.  A  moist  convective 
adjustment  process  is  included  to  represent  cumulus 
convection.  A  staggered  grid  is  used  in  the  horizontal 
because  it  gives  a  better  treatment  of  the  solutions  that 
involve  heating. 

In  Section  A,  the  basic  equations  in  pressure 
coordinates  are  developed  including  the  appropriate  moisture 
terms.  The  convective  adjustment  and  latent  heating  are 
discussed  in  Section  B.  In  Section  C,  the  boundary  layer 


treatment  is  described. 


The  initial  conditions  and 


V4  "V 


tropopause  simulation  are  in  Section  D.  The  details  of  the 
specification  of  model  grids,  grid  arrangements  and  finite 
difference  scheme  are  in  Appendix  A,  B  and  C  respectively. 


A.  BASIC  EQUATIONS 

The  hydrostatic  primitive  equations  in  pressure 
coordinates  may  be  written 


|i  +  v.(*\0  (wV)  +V*  +  fkxV  =  V2V  +  -fp  .  <2'1-1) 


t,.(6V)  =  AeV28  *  4-  fK(p)  |il  +  Ofi  +  0H  +  HC  , 


(2.1.2) 


|iu,.(qv)  ♦  -§-(„,)  =  A„72q  +  4-  CK(p)  |1]  +  M*  -  Mc  ,  (2.1.3) 


V  »V  +  5w/&  p  =  0, 


d<t>^  p  =  -r(p)6  , 


where : 


V  =  horizontal  component  of  the  velocity, 
w  =  dp/dt, 
p  =  pressure, 

<j»  =  gZ  =  geopotential , 

6  =  T(p0/p)K  =  potential  temperature, 

T  =  temperature, 

R  =  gas  constant, 

-p  =  specific  heat  at  a  constant  pressure, 
k  =  R/Cp, 


(2.1.4) 


(2.1.5) 


.  •v-’-.v.v  Z-’Z*  v.-'. 


--  '/  V  >  /  V  .*  -*  ■  *  y 
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f  =  2Q  sin  4-  =  Coriolis  parameter, 

Q  =  rotation  rate  of  earth, 

4>  =  latitude, 
q  =  specific  humidity, 

Am  =  horizontal  diffusion  coefficient  of  momentum, 

Aq  =  horizontal  diffusion  coefficient  of  heat, 

Aq  =  horizontal  diffusion  coefficient  of  humidity, 

K(p)  =  vertical  diffusion  coefficient, 

Qa  =  heating  due  to  dry  convective  adjustment, 

0M  =  heating  due  to  moist  convective  adjustment, 

Hf;  =  large  scale  latent  heating, 

M/\  =  specific  humidity  change  due  to  moist  convective  adjustment, 
Mq  =  specific  humidity  change  due  to  large-scale  condensation, 
r(p,  .  *p-i<-^yc  -  c„ . 


3  5 


The  upper  and  lower  boundary  conditions  are  given  by 


«(x,  y,  Pt»  t)  =u)(x,  y,  p0,  t)  =  0, 

V(x,  y,  p0,  t)  =  0, 

5  V 

eTp(x»  y»  Pt»  ^  =  °»  (2.  i.6) 

|^(xt  y»  Pt.  t)  =  y.  Po*  t)  =  o. 

Here  pt  and  po  represent  the  pressure  at  the  top  and  bottom 
of  the  atmosphere,  respectively,  and  these  quantities  are 
approximated  as  constants.  The  conditions  (2.1.6)  include 
the  kinematic  boundary  condition  at  the  top  and  bottom  of 
the  atmosphere,  the  zero  slip  condition  at  the  surface,  the 
zero  stress  condition  at  the  top  of  the  atmosphere,  and  the 
zero  heat  flux  conditions  at  the  top  and  bottom  of  the 
atmosphere,  respectively. 

The  approximate  steady-state  solutions  to  (2.1.1)- 
(2.1.5)  that  satisfy  the  boundary  conditions  (2.1.6)  are  as 
foil ows : 

ii  =  nrx(l  -  e-^  cos  l )  +  yz<  sin  K 1 

v  =  -nry(l  -  cos  F, )  +  sin 
a)  =0, 

(2.1.7) 

<t>  =  $  =  -D2  (x2  +  y2  )/?  -  fDxy  +F(p), 

9  =  G(p), 

q  =  H(p), 


where  £ =  [  f  /  2K  (  p0  )  ]  1  7  2  (  po  -  p).  The  deformation  D,  and  the 

Coriolis  parameter  f  are  constant,  and  the  functions  Ftp), 
G(p)  and  H(p)  are  arbitrary  except  that  F  and  G  must  satisfy 
the  hydrostatic  equation  (2.1.5).  These  solutions  are  only 
approximate  within  the  Kkman  layer  because  some  small 
advection  terms  have  been  neglected. 

The  dependent  variables  are  subdivided  as  follows: 

'J  =  u(x,  y,  p)  +  u(y,  P,  t) , 

V  =  v(x,  y,  p)  +  v(y,  P,  t),  (2.  1.8) 

w  *«(y,  P,  t),  0=e  (y,  p,  t), 

<t>  =$(x,  y)  +  *(y,  p,  t),  q  =  q(y,  p,  t). 


where  all  of  the  departures  from  the  steady-state  solutions 
are  assumed  to  be  independent  of  x.  If  we  substitute 


the  relations  (2.1.8)  into  (2.1.1) 


2.1.5) 


and  use 


du 

b  X 


d  v 


=  H(1  -  e~^  cos  F, ) ,  =  De-^  sin  u,  v  satisfy 


5  v 


by  ' *  '»  by  5x 

the  steady-state  equations.  The  equations  of  the  departure 

field  from  the  steady-state  solution  become 


djj 

bt 


+ 


(uv)  + 


+  II 


d_u 
5  x 


+ 


aju 

bP 


=  fv  +  A 


b2  u 
m  5  yr 


+  JL.  (x  ?iL) 

bp  5  p; 


* 


(2.1.9) 


a  v  a  ,  ,  ,  a  ,  »  .  a  r  \  j.  a  v  .  a  v 

vt  *  17  ,vv)  *  Tp  luv)  *  17  (vv)  *  “  s—  +  “  5p 


a-ii  .  a2  v  .  0  ov> 

Sy  +  A™  by2'  ~~b$  (K  Sp] 


=  -fu  -  ~T  +  A, 


(2. i . 10) 


ae  ,  a  /  „ »  ,  a  .  a «  ,  -  ae  .  a2e  a  50, 
JT  *  T7  (ve  1  +  T?  (ue '  *  v  57  %  «F  "»7  (K  Sp> 


+  Oa  +  0m  •+  > 


(2.1.11) 


5“t  +  T7  (vq)  +  Tp  (wq)  +  ~v  ly  =  Aq  +  -h  {K  lp} 


+  NA  -  Mc  , 


bs  +  dJ±  =  n 
ay  ap 


(2.1.1: 


1 .  1  :< 


S?  ■ 


.1.14 


Using  the  relations h ip  of  (2.1.7),  the  (2.1.9)  becomes 


5T  +  Ty  (uv)  +  -fp  ^wu)  +  n(1  "  e<  cos  *•  )u 


"  ^  -e-^  cos  F, )  -  xe"^  sin  ^  1  —  +  D(e"^  sin  Z  )v  (2.1.15 

+  wn  [x(l  -  cos  Z)  +  ye<  sin  ?  ] 

=  fv  +  Am  +  -I—  (K  — ) 

m  dy2  ap  v  ap;  ' 


s.v.V. 


.  V’  -\.  \-\.W'  V  ■V-  *\. 


It  is  to  be  noticed  that  the  assumption  that  u  is 
independent  of  x  for  all  time  is  violated  in  (2.1.15)  in  the 
two  terms  where  x  appears.  This  indicates  that  the  other 
quantities  which  were  assumed  independent  of  x  are  actually 
x-dependent,  which  would  bring  other  terms  into  the 
equations.  In  this  study,  these  added  terms  are  neglected 
and  (2.1.15)  and  the  other  equations  are  applied  at  x=0.  If 
the  initial  u,  v,  a),  0  and  q  fields  are  independent  of  x, 
the  error  from  this  approximation  will  grow  slowly  and  it 
will  be  confined  to  the  boundary  layer  since  the  x-dependent 
terms  in  (2.1.15)  are  zero  outside  the  boundary  layer. 
Since  frontogenesis  occurs  very  rapidly,  it  is  expected  that 
the  development  of  x-variations  in  the  dependent  variables 
would  only  have  a  small  effect  on  the  resulting  quasi-steady 
front.  In  any  case,  these  effects  could  not  be  observed  in 
the  atmosphere,  because  atmospheric  fronts  always  have  some 
variation  in  tiie  basic  fields  along  the  front  which  would  be 
much  more  important. 

When  (2.1.15)  is  evaluated  at  x  =  0,  it  becomes 

U  +  +V"U)  +n(1  -e^  cns^+r!y 

+  n(e^  sin  I ) v  -  uiVn  sin  5  )  =  fv  +  Am  ~r  +  ~  (K  ~) 


(2.1.16) 


V  r‘  Ftf 


where  VD  =  -Hy  and  T  =  Vn(l  -  e"S  COS  l).  To  remove  the 

vertical  mean  in  (2.1.10),  the  vertical  average  is  defined 


<(  )>  = 


Po  "  Pt 


Po 

/  (  )hp. 

Pt 


By  vertically  averaging  (2.1.13)  and  using  the  boundary 
conditions  (2.1.6),  we  obtain 


<v>  =  n  . 


This  equation  states  that  the  total  disturbance  y  mass  flux 
is  independent  of  y.  It  can  be  expected  that  at  a  large 
distance  from  the  frontal  gone  the  disturbance  mass  flux 


will  vanish.  Thus,  we  get 


<V>  =  0  . 


(2.1.17) 


If  we  take  the  vertical  average  of  (2.1.10)  and  use  (2.1.6) 
and  (2.1.7),  we  obtain 


ly  <W>  +  <"W>  *  <»  +  ^  ^ 

,  .  .  f<u>  + 


(2.1.18) 


\T 


'-M  •*;  V,' fj  v; v,-  v;v; g-;y  ,■  V. ■g-j’j  .  ■-.  T-.  ir.  -y  ■ 


When  this  result  is  subtracted  from  ( 2*.  1.10),  we  have 


U  *  T7  (vv  •  <vv>)  +  Tp  (“v)  +  T7  (~vv  ■  <''v>) 


a v  .  a v.  ,  av  .  a  v. 
*  "  57  -  <u  St  *  “  dp  ‘  *  S? 


=  -  f(u  -  <v>)  — (it  -  <rt>) 

♦  V  Sjf ♦  -fc  »U)  ■  (4-'B 


ap  v  a p * 


p0  SP; 


Po 


Using  (2.1.7)  and  evaluating  at  x  =  0,  we  obtain 


It  *  T7  («  -  <»*»  +  -§7  <“V>  *  T7  <rv  -  <rv» 


+  n[ue“£  sin  £  -  <ue"£  sin  £>  -  V[)[oj  (e*^  cos  'i ) 
9  (n  -  <n>)  -  f(u  -  <v>)  +  A„, 


{e<  cos  &)>] 


ay 


+  JL_  (K  ll)  .  (JL  11) 

ap  (K  ap}  Po  Qp 


(2.1.19) 


Po 


Eq. (2.1.11)  evaluated  at  x  =  0  may  be  written, 


ae  x  a(ve)  a  ,  *  ae 
at  ^y”  +  Tp  (we)  +r  5J 


(2.1.20) 


=  Ae  +  —■  TK(p)  ||]  +  nA  +  0M  +  HC 


A  ^ 


Similarly,  (2.1.12)  becomes 


£t  +  Ty  (vp)  +  Tp  +r  Ey 

2  (2.1.21) 

=  Aq  ay^  +  Tp  tK(p)  !p]  +  ma  -  MC  • 

Integrating  the  hydrostatic  equation  (2.1.14)  with 
respect  to  p  and  removing  the  vertical  mean,  we  obtain 

,P  P 

tc  -  <rr>  =  J  r(p)9dp  -  <  /  r(p)0dp  >  .  (2.1.22) 

Pt  Pt 

With  (2.1.13),  (2.1.16)  and  ( 2 . 1 . 19 ) - ( 2 . 1 . 22 )  form  a 

complete  set  of  equations  which  can  be  solved  by  a  marching 
process. 

R.  CONVECTIVE  ADJUSTMENT  PROCESSES  AND  LATENT  HEATING 

This  section  describes  the  dry  and  moist  convective 
adjustment  processes  and  the  latent  heat  release  that  are 
used  in  the  model.  The  dry  convective  adjustment  process  is 
activated  in  any  layer  where  99  /  9p  >  0.  In  such  a  layer, 

0  is  set  equal  to  the  average  value  within  the  layer.  Then 
the  layers  above  and  below  are  tested  for  instability  and 
modified  if  necessary.  This  process  is  represented  by  Q a  in 
(2.1.20). 

When  the  relative  humidity  exceeds  100%,  that  is 
q  >  qs  [ qs ( T , p )  is  the  saturation  mixing  ratio],  latent  heat 


a 
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of  condensation  will  be  released.  To  adjust  the  relative 


humidity  to  100%,  T  and  q  must  be  changed  to  satisfy 


q  +  6q  =  qs(T  +  6T,  p) , 


Cp6T  =  -L6  q , 


(2.2.  1) 


(2.2.2) 


where  L  is  the  latent  heat  of  condensation,  and  6  q  and  ST 
are  the  changes  in  q  and  T  during  the  process.  A  good 
approximation  to  the  function  qs(T,p)  is  given  by 


%  = 


exp  [-5—  (- 


(2.2.3) 


where  Rv  is  the  gas  constant  for  water  vapor. 

The  equations  (2.2.1)  and  (2.2.2)  require  that  the  final 
state  of  the  air  be  exactly  saturated  and  that  excess 
moisture  be  condensed  isobarically,  which  releases  the 
latent  heat  to  the  air. 

Equations  (2.2.1)  and  (2.2.2)  can  be  expressed  in 
iterative  form  as 


6q  -  qs  (1  -  r)  +  (  ST.  H  =  rcls 


CpST  =  -L6q  . 


Eliminating  ST  from  (2.2.3)  and  (2.2.4)  gives 


(2.2.4) 


(2.2.5) 


^ K.‘ V  V  «'  V  v  .*  *.*  *«■  V  *.■  ,*  *,•  ',* 


6q 


q<u  -  r) 


1  + 


_L_  (_La s_) 

-  '  aT  ; 


(2.2.6) 


At  constant  pressure,  the  C laus i us-C lapey ron  equation  can  be 
expressed  as 


1  dT 


■) 


Substituting  (2.2.6)  into  (2.2.7)  gives 


(2.2.7) 


qg(l  -  r) 

6fl  =  i2n 

1  +  "c^- 

S  q  may  now  be  calculated  from  (2.2.8)  as  can  the 
corresponding  value  of  6T  from  (2.2.4).  The  first  estimates 
of  the  adjusted  temperature  and  mixing  ratio  are, 

T'=T+6T  and  q'=q+6q.  (2.2.9) 

These  values  are  in  turn  used  to  obtain  improved 

approximations.  Iterations  are  continued  until  the  values 
of  T’  and  q’  at  each  level  match  with  sufficient  accuracy  a 
value  of  100%  relative  humidity.  The  changes  are  then 

applied  to  the  potential  temperature  and  mixing  ratio  fields 


of  the  model,  which  accounts  for  the  term  He  of  (2.1.20) 


The  convective  latent  heating  is  parameterized  by  a 


convective  adjustment  scheme.  The  scheme  yields  similar 
results  in  describing  the  large-scale  effect  of  latent 
heating  as  an  explicit  scheme  which  included  prediction 
equations  for  water  vapor,  cloud  water  and  rain  water  (Hsie 
and  Anthes,  1904). 

In  areas  where  convection  is  occurring,  it  is  necessary 
to  determine  whether  the  atmosphere  is  convectively  unstable 
to  moist  processes.  A  layer  will  be  conditionally  unstable 
if  39e / 3p  >  0,  where  0e  is  the  equivalent  potential 
temperature,  which  is  given  by 

0e  =G  exp(— )  .  (2.2.10) 

tp  i 

As  with  the  dry  convective  adjustment,  0e  is  set  equal  to 
the  average  value  in  the  unstable  layer.  This  gives  a  value 
for  0p  at  each  level  and  6T  and  <5  q  can  then  be 
determined  from  the  following  calculations. 

The  heat  release,  Qm  ,  as  a  result  of  moist  convective 
adjustment  is  included  through  knowledge  of  the  change  of 
equivalent  potential  temperature  during  the  adjustment. 
Taking  the  natural  logarithm  of  (2.2.10)  gives 

Jin  8e  =  in  9  + 

Cp  i 


4  5 


(2.2.11) 


A  constant  pressure  surface  may  be  assumed  since  the  model 


uses  pressure  coordinates.  Differentiating  both  sides  of 
(2.2.11)  and  recalling  the  definition  of  potential 
temperature  gives 


60p  6  T  .  6 dcL 

6e  "  T  CpT 


(2.2. 12) 


Solving  (2.2.6)  for  dqs ,  substituting  into  equation  (2.2.11) 
and  solving  for  dT  leads  to 


6  T  = 


60, 


0e  yCU  + 


q*1-2- 


RyCpT 


J=3s_ 

cpT 


(2.2.13) 


where  50  e  is  the  change  in  0e  needed  to  adjust  the 
gradient  30e/3p  to  be  equal  to  zero. 

Since  the  atmosphere  remains  saturated  throughout  the 
process  of  moist  convective  adjustment,  the  Clausius- 
Clapeyron  equation  is  used  again  to  obtain  the  corresponding 
6  q  or  actual  5qs.  Changes  are  again  applied  to  the 
potential  temperature  and  mixing  ratio  fields  of  the  model 
which  accounts  for  the  terras  Qm  of  equation  (2.1.20)  and  Me 
of  equation  (2.1.21).  The  convective  adjustment  with  a 
saturated  atmosphere  is  a  rather  crude  parameterization  of 
the  cumulus  effects.  However,  the  relative  humidity  can 
indeed  become  very  high  in  intense  Mei-Yu  cases  (Chen  and 
Tsay,  1978).  Furthermore,  the  strong  control  by  large-scale 


§ 
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f rontogenes is  forcing  makes  the  results  less  sensitive  to 
the  details  of  the  cumulus  parameterization. 


C.  BOUNDARY  LAYER  PARAMETERIZATION 

The  surface  fluxes  are  parameterized  by  the  bulk- 
aerodynamic  formula  and  the  surface  fluxes  of  heat  and 
momentum  are  assumed  to  influence  directly  only  the  lowest 
model  layer.  The  magnitude  of  the  surface  stress  is 
proportional  to  the  square  of  the  wind  speed  and  the  surface 
sensible  heat  flux  is  proportional  to  the  product  of  the 
wind  speed  and  the  vertical  gradient  of  temperature  between 
the  ground  and  the  middle  of  the  lowest  model  layer.  This 
formulation  is  best  suited  to  models  with  coarse  vertical 
resolution,  since  the  depth  of  the  lowest  layer  should 
correspond  to  a  typical  PBL  depth. 

In  this  study,  the  vertical  sub-grid  scale  mixing  terms 
in  Eqs. (2.2.16)  are  applied  to  the  lowest  layer  of  the 
model.  The  term  may  be  written  3s 


-Spowlfi 


At  the  top  of  the  boundary  layer,  we  assume 


(2.3.1) 


F„(p)  =  =  K<z) 


a  i  * 


(2.3.2) 


K(z)  is  defined  as 


K(x)  -  cD j u | 


where  Cd  is  the  drag  coefficient  and  |u|  =  t/u2  +  V2  . 
Applying  (2.3.2)  to  the  lowest  two  layers,  po  and  pi ,  of  the 


model,  we  get 


Fu  =  -  pg  Cn|uj(ui  -  uq)/(pi  -  po) 


=  -  pg  Cn ) u | ux /A p 


since  uo  =  0  from  the  lower  boundary  condition  and  Ap  =  p-j  - 
Therefore,  the  u  momentum  flux  at  the  surface  becomes 


fu  =  -  in-3  ( 


•)  ui  * 


(2.3.3) 


where  10~  3  is  a  unit  conversion  factor  which  converts  the 
unit  of  pressure  to  cb .  Similarly,  the  v  momentum  flux  is 


F«  =  -  10-3 


(2.3.4 


In  this  study,  the  PBL  depth,  h,  is  estimated  as 


Pn  -  Pi 


mo  -91.9  in3 

1.25  x  9.81  x  10  =  66(1  m* 


The  sensible  heat  and  moisture  fluxes  are 


pe  -  -  in-3  (e  x  -  »0) 


(2.3.5) 


Fq  .  .  m-3  (qi  -  qn)i 


(2.3.6 
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where  9o  and  qo  are  the  specified  surface  values 
The  latent  heat  flux  is  defined  as 


Ft  - 


L  »> 


(2.3.7) 


where  L  =  2.5  x  10s  Joule/kg,  cp  =  1003  Joule/ (° K/kg ) . 

D .  INITIAL  CONDITIONS 

The  initial  potential  temperature  is  given  by 


2 

0(y,  P»  n)  =  e  I  ( p )  -  a(-^-)  arctan  (sinh  ay)  , 


(2.4.1) 


where  ©  1  ( p )  is  a  prescribed  vertical  potential  temperature 
profile  and  a  is  one-half  the  total  horizontal  potential 
temperature  variation.  a  is  defined  as  and  is  evaluated 


frt/H  v-~—  (— ( -  ~v~M  at  some  middle  level  p  =  pm  . 

Pm  Pm  3  p  Pm 

The  thermal  wind  equation,  which  is  derived  using  the 
hydrostatic  equation  and  the  definition  of  geostrophic  wind, 
may  be  written  as 


-  II  (_£_)<  ae_ 
aP  f  p  1  Po  ay  * 


(2.4.2) 


Integrating  (2.4.2)  from  p  =  p™  to  p  =  p 


rp  Ly.  dn  -  R  50  rP  1  P  ^  ^ 

1  »P  dp  TSy  p/  p  “P 


Pm 


and  assuming  u  =  0  at  p  -  pm ,  we  obtain 
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p..*- 


u(y,  P.  0)  -  ^  C(J_f  -  (Mfl  . 

Using  (2.4.1), 

57“  ■  ^TT7  rarcUn  (Sinh  ay)7 
=  -  sech  (ay)  ■. 

it 


(2.4.3) 


(2.4.4) 


Substituting  (2.4.4)  into  (2.4.3),  the  initial  u  component 


of  the  velocity  may  be  written, 


u(y,  p,  0)  =  ~  ~ ((-£-)*  -  sech  (ay)  .  (2.4.5) 

'  T*-  Pn  Dr» 


The  v  component  of  the  velocity  is  set  initially  to  zero, 
i  .  e.  , 

v(y,  p,  0)  =  0  , 

so  that ,  w  (y,  p,  D)  =  0  . 

The  initial  field  of  mixing  ratio  is  calculated  by 
determining  the  saturation  vapor  pressure  es  at  each  grid 
point  from  an  integrated  form  of  the  C 1 aus i us-C 1 apey r on 
equation.  According  to  Hess  (1959), 

es  =  0.611  exp  (-Jjr-  (~y^  -  j  ))  »  (2.4.6) 


where  To 

is  a 

constant 

of  integration. 

Knowing  the 

siiturat  ion 

v  a  p  o  r 

p  i  <  •  s  s  u  v  e  , 

the  so  t  u  rn  t  i  on 

in  i  x  i  n  g  ratio  i  s 

E 


obtained  from  the  approximate  relation  (Haitiner  and  Martin, 


1957)  . 


0.622  ec 

qs  p 


(2.4.7) 


Knowing  the  potential  temperature  at  t  =  0  from  (2.4.1),  it 
is  possible  to  calculate  qs  at  each  point  using  equations 
(2.4.6)  and  (2.4.7).  Defining  q  to  be  a  given  fraction  of 
qs  allows  the  testing  of  atmospheres  with  different  initial 
relative  humidities. 

The  initial  conditions  described  in  previous  sections 
will  be  constrained  by  a  simulated  tropopause  that  is 
centered  at  about  300  mb  level.  The  simulation  is 
accomplished  by  using  a  reference  temperature,  T ref,  at 
about  the  midpoint  of  the  horizontal  scale  and  calculating 
temperatures  such  that  a  tropopause  with  the  expected 
variation  in  height  is  determined  and  a  zero  vertical 
temperature  gradient  is  achieved  in  the  lower  stratosphere. 
A  reference  temperature,  Tref,  is  computed  using  (2.4.1) 


Tref  =  215  +  a 


)  arctan  (sinh  ay) 


(2.4.8) 


When  the  temperature  in  each  column  reaches  this  reference 
temperature,  it  is  held  constant  at  that  value,  as  expected 
in  the  lower  stratosphere.  The  temperature  is  then 


converted  to  potent!  a  l  temperature. 


III.  NUMERICAL  EXPERIMENTS  AND  RESULTS 


In  this  chapter  two  groups  of  numerical  experiments  are 
carried  out  to  simulate  the  midlatitude  and  subtropical 
situations,  respectively.  For  each  group,  the  Coriolis 
parameter  and  initial  temperature  sounding  are  fixed,  but 
the  amount  of  the  initial  relative  humidity  and  surface 
fluxes  of  momentum,  mixing  ratio  and  sensible  heat  are 
varied.  The  purpose  here  is  to  describe  the  structure 
differences  due  to  the  difference  in  latitude,  relative 
humidity  and  surface  fluxes.  The  discussion  of  the  results, 
which  focuses  on  the  development  of  the  low  level  jet,  is 
giveninChapterIV. 

We  will  first  give  the  parameter  ranges  used  in  the 
experiments  (in  Section  A),  and  then  describe  the  numerical 
results  of  the  midlatitude  (Section  B)  and  subtropical 
(Section  C)  cases,  respectively.  As  mentioned  in  the 
Introduction,  the  main  purpose  for  all  numerical  experiments 
is  to  obtain  a  quasi-steady  solution.  However,  it  is  also  of 
interest  to  follow  the  time  evolution  of  the  frontogenesis 
process  to  see  how  the  quasi-steady  states  are  developed. 
Therefore,  in  each  group  the  evolution  for  a  "typical  case" 
is  shown.  It  turns  out  that  these  evolutions  provide  useful 
insigiits  about  the  different  nature  of  the  effect  of  cumulus 
convection  in  the  midlatitude  and  subtropical  cases. 


A.  EXPERIMENT  CASES  AND  PARAMETERS 

The  following'  values  are  used  for  the  constants  in  this 

s  t  udy : 


=  9.81  m  s"2  . 

=  1  Of) 3  joul e/kg/°K 
=  287  jou1es/kg/°K 
=  R/Cp 

=  461  joules/kg/°K 
=  2.5  x  lf£  joules/kg 

=  5°K/km 


(^)pm  =  “  0.53448°K/cb 
p0  =  100  cb 


Pm 

=  55  cb 

Am 

=  Aq  =  ^ 

=  3  x  104  n?/s 

s? 

II 

cr 

ii 

=  0.0007  cb2/s 

At. 

=  270  sec 

Ay 

=  60  km 

Y 

=  2760  km 

Each  experiment  is  designated  by  a  three-letter  case 
name  as  follows: 


LRF 

where  L  defines  the  latitude  environment. 

1.  L  =  M  refers  to  a  midlatitude  environment,  with  the 

following  parameter  settings: 

a.  f  =  1.06  x  10-4  s-1  (corresponds  to  46°N); 

b.  initial  vertical  potential  temperature 

distribution  given  by  a  typical  midlatitude 

sounding  (Fig.  3.1); 

c.  initial  meridional  temperature  variation  is 

given  by  a  =  8°K;  and 

d.  large-scale  deformation  fields  given  by  D  = 

1  x  10-s  s'1  . 

2.  L  =  S  refers  to  a  subtropical  environment,  with  the 

following  parameter  settings: 

a.  f  =  0.558  x  10-4  s"1  (corresponds  to  22.5°N); 

b.  initial  vertical  potential  temperature 

distribution  given  by  a  typical  subtropical 

sounding  ( Fig .  3.1); 

c.  initial  meridional  temperature  variation  is 

given  by  a  =  4°K;  and 

d.  large-scale  deformation  fields  given  by  D  = 


0.5  x  10'5  s'1. 


3.  R  defines  the  relative  humidity: 
a.  R  =  6  means  R.H. =60%;  and 


b.  R  =  9  means  R  .  H  .  =  9  0  %  . 

The  observed  low-level  relative  humidi'y  in  the  vicinity 
of  a  Mei-Yu  front  varies  between  60%  -  90%  (Chen  and  Tsay, 
1978).  When  the  Mei-Yu  front  is  north  of  the  Yantze  River  (a 
midlatitude  situation)  the  minimum  R.H.  is  about  60%.  When 
it  is  located  over  Taiwan  and  southeastern  coast  of  the 
Chinese  mainland  the  maximum  R.H.  is  near  90%. 

4.  F  defines  the  surface  fluxes: 

a.  F  =  0  means  there  is  no  surface  fluxes  of  heat, 
moisture,  and  momentum;  and 

b.  F  =  1  means  surface  fluxes  of  heat,  moisture  and 
momentum  are  included. 

Typically,  the  fluxes  are  large  when  a  subtropical 
portion  of  the  front  is  over  a  warm  sea  surface.  When  the 
front  is  moved  over  land  in  the  midlatitude  cases,  the 
fluxes  are  normally  weaker  with  the  exception  of  momentum 
fluxes.  For  self-consistency,  the  fluxes  of  all  quantities 
will  either  be  included  or  excluded  in  the  numerical 
experiments.  For  all  cases,  the  initial  zonal  wind  is 
calculated  from  the  potential  temperature  distribution  using 
the  thermal  wind,  and  the  initial  v  component  of  the 
velocity  is  set  to  zero. 


development  of  the  frontal  thermal  structure,  including  the 
concentration  of  the  isentropes  in  the  lowest  2  km  and  the 
poleward  tilt  of  the  frontal  zone  from  the  surface  to  the 
upper  troposphere,  is  clearly  indicated  in  Fig.  3.4.  In 
addition,  a  secondary  maximum  occurs  in  the  upper 
troposphere  (6  -  8  km  or  600  -  400  mb)  near  50°,  and  a 
reversed  tilt  and  reversed  gradient  is  present  in  the  lower 
stratosphere.  This  stratospheric  structure  is  a  result  of 
the  imposed  condition  of  a  tropopause  and  is  consistent  with 
previous  results  obtained  by  Williams  et  al.  (1981). 

The  vertical  velocity  field  (Fig.  3.5)  quickly 
develops  and  becomes  organized,  with  a  concentrated  upward 
motion  of  >  -4.0  x  10~4  cb/sec  occurring  at  day  2  in  the 
lower  troposphere  near  43°.  A  poleward  vertical  tilt  of  the 
area  of  upward  motion  develops  and  becomes  more  obvious  with 
time.  After  day  4,  a  quasi-steady  state  clearly  delineates 
both  the  upright  maximum  rising  motion  centered  at  43°  to 
44°  in  the  lower-middle  troposphere  and  the  vertical  tilt, 
which  is  well  represented  by  the  0)  =0  line  poleward  of  the 
maximum  rising  motion. 

A  very  small  area  with  a  secondary  maximum  in 
vertical  velocity  appears  at  day  2  near  the  surface  ahead  of 
the  front  (41°).  This  represents  the  effect  of  the  forcing 
by  the  Ekman  layer  structure  of  the  mean  flow,  as  was 
discussed  by  Williams  et  al.  (1981).  Because  this  area  is 
very  close  to  the  surface,  it  does  not  cause  condensation 


and  latent  heat  release  even  in  the  highly  moist  cases 
discussed  below. 

The  initial  zonal  wind  (Fig.  3.2b)  is  in  geostrophic 
balance  with  the  temperature  field,  with  an  upper 
tropospheric  westerly  maximum  of  ~55  m/s  and  a  surface 
easterly  maximum  of  -25  m./s.  The  pattern  is  almost 
symmetric  with  respect  to  46°  latitude.  As  the  integration 
proceeds  (Fig.  3.6),  asymmetries  develop  and  a  strong  shear 
begins  to  appear  at  day  1.0  below  6  km.  The  shear  quickly 
reaches  a  steady  state  with  a  frontal  structure  which  tilts 
from  42°  at  the  surface  towards  45°  at  4  km.  This  is  in  the 
same  vicinity  as  the  tu  -0  line  and  the  frontal  zone  that  is 
defined  by  the  potential  temperature  and  mixing  ratio 
gradients.  At  the  latitude  of  the  surface  front,  a  local 
westerly  maximum  develops  near  z  =  4  km.  By  day  2  (Fig. 
3.6b),  this  narrow  westerly  regime  penetrates  downward  and 
appears  to  be  an  intrusion  into  the  lower  tropospheric 
easterlies.  It  remains  in  the  3-4  km  layer  through  day  4. 
Although  the  magnitude  of  westerly  wind  in  this  narrow  zone 
is  small,  the  poleward  meridional  shear  is  large,  which 
results  in  a  vorticity  in  this  region  that  is  comparable  to 
the  values  near  the  upper  tropospheric  westerly  jet.  The 
development  of  this  narrow  westerly  regime  is  a  result  of 
the  intensifying  meridional  temperature  gradient  which  is 
forced  by  the  f  r  on  t  ogenes  i  s  process  in  the  lower 
troposphere.  Development  of  this  zonal  velocity  field  is 


also  dramatically  illustrated  by  the  evolution  of  the 
relative  vorticity  (Fig.  3.7).  The  initial  flow  pattern 
develops  rapidly  into  a  quasi-steady  state  frontal  structure 
after  about  two  days.  Maximum  cyclonic  vorticity  in  the 
lower  troposphere  is  concentrated  on  the  warm  side  of  the 
front.  To  illustrate  the  vertical  structure  of  the 
horizontal  wind  shear  with  that  of  the  horizontal 
temperature  gradient,  3 0 /  9y  =  -2°K/100  km  for  day  4  is 
plotted  in  the  vorticity  figure.  A  coherent  vertical  tilt 
of  the  frontal  structure  is  obvious. 

2.  Other  Midlatitude  Cases 

The  effects  of  varying  the  initial  relative  humidity 
(M90),  and  the  inclusion  of  surface  momentum,  heat  and 
moisture  fluxes  (M61  and  M91)  are  illustrated  by  the 
resulting  structure  at  day  4  (Figs.  3.8-3.12).  In  addition, 
the  precipitation  and  heating  rates  are  also  shown  in  Figs. 
3.13  and  3.14,  respectively. 

a.  Effect  of  Moisture 

(1)  No  Surface  Flux  Cases  (M60  vs.  M90) 

It  is  clear  that  the  initial  R.H.=90%  case 
(Fig.  3.8b)  has  higher  mixing  ratios  in  general  and  a  warmer 
upper  troposphere  between  40°-50°  (best  represented  by 
comparing  the  0=35O°K  line).  The  higher  temperature  is  a 
result  of  a  larger  latent  heat  release  in  the  frontal  region 
due  to  the  increased  moisture  content,  as  can  be  seen  by 
comparing  the  precipitation  rates  (Fig.  3.13).  The 


meridional  potential  temperature  gradient  in  the  surface 
front  for  the  R.H.=90%  case  is  also  stronger  (Figs.  3.9a 
versus  3.9b).  The  frontal  zones  as  defined  by  the  region 
38/3y  <  -0.02°K  km-1  are  quite  similiar  in  width  and 
vertical  tilt.  Above  the  surface  front  position  (~42°),  the 
mixing  ratio  in  the  R.H.=90%  case  (Fig.  3.8b)  also  shows  a 
significant  small  scale  structure  with  a  wavelength  near  2  y 
in  the  2-7  km  layer.  This  structure  is  much  less  evident  in 
the  R  .  H .  =  6  0  %  case  (Fig.  3.8a).  This  small  scale  structure 
is  near  the  region  of  maximum  rising  motion  above  the 
surface  front  position.  The  stronger  vertical  velocity  for 
the  R.H.=90%  case  (maximum  rising  motion  >  6x10  4  cb/s 
versus  4xl0-4  cb/s  for  the  R.H.=60%  case)  is  consistent  with 
the  stronger  cumulus  convection.  This  convection  arises 
because  of  the  existence  of  conditional  instability  which  is 
due  to  the  prescribed  90%  relative  humidity.  The  cumulus 
convection  activity  appears  to  cause  the  small-scale 
variations  in  the  mixing  ratio  field. 

In  addition  to  a  larger  maximum,  the 
vertical  velocity  of  the  R . H . =  9  0  %  (Fig.  3.10b)  case  also 
displays  a  more  complex  structure,  with  two  concentrated 
regions  of  rising  motion.  The  primary  region  near  42°-43° 
extends  throughout  most  of  the  depth  of  the  troposphere  from 
about  1  km  to  11  km.  The  secondary  region  occurs  in  the 
middle-upper  troposphere  (5-11  km)  just  poleward  of  the 
primary  region  at  about  44°.  In  addition,  a  third, 

60 


relatively  weak,  maximum  upward  motion  occurs  farther 
poleward  near  47°.  This  structure  is  referred  to  as  a 
multi-band,  slanted  rising  motion  in  later  discussions.  This 
multi-band  structure  in  the  vertical  velocity  may  be 
compared  to  some  of  the  observed  multi-band  structure  of 
cloud  and  precipitation  associated  with  cold  fronts  during 
the  period  of  CYCLES  (CYCLonic  Extratropical  Storms)  Project 
(Parsons  and  Hobbs,  1983). 

As  may  be  expected,  the  stronger  warm  core 
and  stronger  convection  of  the  R.H.=90%  case  is  related  to  a 
larger  latent  heat  release.  In  Fig.  3.14,  the  heating 
distribution  generally,  resembles  the  rising  motion  field 
(Fig.  3.10).  This  is  consistent  with  a  near  balance  between 
adiabatic  cooling  and  diabatic  heating.  The  R.H.=90%  case 
(Fig.  3.14b)  has  a  maximum  heating  rate  of  >6  °K/day,  which 
is  much  higher  than  the  1  °K/day  rate  of  the  R.H.=60%  case 
(Fig.  3.14a).  This  six-fold  increase  in  heating  is 
consistent  with  the  increase  in  precipitation  rate  shown  in 
Fig.  3.13.  The  greater  production  of  kinetic  energy  implied 
by  the  stronger  rising  motion  of  warm  air  gives  rise  to  a 
more  intense  horizontal  circulation  (Figs.  3.11b  versus 
3.11a).  The  upper  tropospheric  westerly  jet  (Fig.  3.11b) 
reaches  a  maximum  of  >  50  m/s  (versus  the  >  35  m/s  for  the 
R.H.=60%  case)  near  52° ,  which  is  two  degrees  poleward  of 
the  jet  core  of  the  R  .  H .  =  6  0  %  case.  A  small  region  of  weak 
easterlies  that  develops  equator ward  of  the  westerly  jet  is 


not  observed  in  the  R.H.=60%  case.  It  is  to  be  recalled 
here  that  these  are  perturbation  wind  fields;  the  basic  flow 
is  westerly  throughout  the  upper  troposphere.  The 
distribution  of  the  perturbation  zonal  wind  is  consistent 
with  an  anticyclonic  circulation  above  the  warm  core  that  is 
due  to  heating.  The  existence  of  this  easterly  regime 
limits  the  middle-upper  tropospheric  westerlies  to  poleward 
of  the  surface  front  position,  and  makes  the  middle 
tropospheric  westerly  regime  appear  as  a  narrow  jet-like 
structure  that  penetrates  downward  to  below  z-3  km.  The 
resultant  low  level  westerly  regime  is  sharper  and  more 
extensive  than  that  of  the  R.H.=60%  case. 

On  the  other  hand,  the  strength  of  the 
surface  easterlies  of  the  M60  and  M90  cases  is  very  similar. 
This  is  probably  because  most  of  the  impact  of  a  higher 
relative  humidity  is  in  deep  cumulus  convection,  with  the 
resultant  increased  latent  heat  release  occurring  mainly  in 
the  upper  and  middle  troposphere,  as  shown  in  Fig.  3.14. 

The  maximum  cyclonic  vorticity  is 
concentrated  within  the  frontal  zone  in  both  the  upper  and 
the  lower  troposphere  for  both  cases  (Fig.  3.12a  versus  Fig. 
3.12b).  The  stronger  horizontal  circulation  in  the  upper 
troposphere  for  the  R.H.=90%  case  is  also  evident  in  the 


vorticity  field. 


(2)  Surface  Flux  Cases  (M61  versus  M91) 

When  surface  fluxes  of  momentum,  moisture 
and  heat  are  present,  the  differences  between  the  more  moist 
(R.H.=90%)  and  the  drier  (R.H.=60%)  cases  are  basically 
similar  to  those  between  the  no-flux  cases.  The  R.H.=90?o 
case  has  a  warmer  upper  tropospheric  core,  as  represented  by 
the  decrease  of  height  of  the  0  =350°K  line,  and  again  has 
the  small-scale  variation  in  mixing  ratio  in  the  rising 
motion  region  over  the  surface  front  location  (Fig.  3.8d 
versus  Fig.  3.8c).  The  meridional  potential  temperature 
gradient  is  stronger  for  case  M91  than  for  case  M61  (Fig. 
3.9d  versus  Fig.  3.9c).  The  maximum  upward  motion, 
precipitation  and  heating  rate  are  much  larger  for  M91  than 
for  M61  (Figs.  3.10,  3.13  and  3.14).  Stronger  wind 
circulations  (Fig.  3. lid  versus  Fig.  3.11c)  are  shown  in  the 
R.H.=90%  case.  The  circulation  includes  a  more  intense 
westerly  regime  of  >  50  m  s_I  in  the  upper  troposphere,  the 
development  of  a  narrow  easterly  regime  equatorward  of  it, 
and  a  deeper  downward  extension  of  the  upper  westerlies  into 
the  low-level  easterly  jet.  The  differences  in  the 
magnitudes  of  vorticity  (Fig.  3.12d  versus  Fig.  3.12c)  are 
similar  to  other  fields  with  much  stronger  cyclonic 
vorticity  within  frontal  zone  for  the  R  .  H  .  =  9  0  %  case. 

It  is  to  be  noted  that  the  small-scale 
variations  in  all  fields  for  the  case  with  high  mixing  ratio 
and  surface  fluxes  (Fig.  3.8d-3.12d)  exist  near  surface  to 


the  north  of  front.  As  can  be  seen  in  Fig.  3.15,  this  is 
the  cold  air  region  where  the  sensible  and  latent  heat 
fluxes  are  large  and  apparently  cause  a  shallow  layer  of 
conditional  instability. 

On  the  other  hand,  there  is  one  noticeable 
difference  between  the  comparisons  of  the  with-flux  cases 
and  those  of  the  no-flux  cases.  The  multi-band  vertical 
velocity  structure  in  the  no-flux  (R.H.=90%)  case  (Fig. 
3.10b)  is  not  as  conspicuous  in  the  with-flux  (R.H.=9Q%) 
case  (Fig.  3.10d).  However,  a  weak  secondary  maximum  region 
does  appear  in  the  latter  near  7  km  at  46°,  which  indicates 
a  double-maximum  structure  and  a  poleward  tilt  in  the 
vertical  structure.  The  physical  implication  of  this 
slantwise  structure  in  the  upward  motion  field  will  be 
discussed  later  in  Chapter  4. 

b.  Effect  of  Surface  Fluxes 

(1)  R.H.=60%  Cases  (M60  vs.  M61) 

The  inclusion  of  surface  fluxes  (Figs.  3.8a 
and  3.8c)  does  not  lead  to  any  discernible  differences  in 
potential  temperature  and  mixing  ratio  distributions,  except 
the  mixing  ratio  isolines  at  the  surface  have  been  spread 
towards  cold  air  when  surface  fluxes  are  included. 
Comparison  between  Figs.  3.9a  and  3.9c  shows  that  the 
surface  meridional  potential  temperature  gradient  for  the 
case  with  fluxes  substantially  weakens,  although  the  general 


pattern  delineating  the  frontal  zone  remains  the  same. 


These  reductions  in  the  potential  temperature  and  mixing' 
ratio  gradients  at  the  surface  are  due  to  the  distribution 
of  the  surface  fluxes  given  by  the  bulk  formula,  which 
provides  larger  heat  and  moisture  fluxes  poleward  of  the 
surface  front  in  the  cooler  and  drier  region  (Fig.  3.15). 

The  inclusion  of  surface  fluxes  causes  a 
slight  increase  in  the  vertical  velocity  (Fig.  3.10a  and 
3.10c).  This  is  expected  since  more  latent  heat  release 
occurs  with  the  increased  supply  of  heat  and  moisture  from 
below,  as  shown  in  the  comparison  between  Fig.  3.14c  and 
Fig.  3.14a.  Since  the  resultant  upper  tropospheric 
temperature  is  not  higher  (Figs.  3.8a  and  3.8c),  the 
additional  latent  heat  release  in  the  area  of  rising  motion 
is  compensated  almost  entirely  by  the  increased  adiabatic 
cooling  associated  with  the  stronger  upward  motion. 

Increased  conversion  of  available  potential 
energy,  as  implied  by  the  in-phase  relationship  between  the 
vertical  velocity  and  potential  temperature  fields,  leads  to 
a  stronger  upper  tropospheric  jet  (40  m/s  versus  35  m/s  in 
the  no-flux  case)  (Figs.  3.11a  and  3.11c).  On  the  other 
hand,  a  decrease  of  the  surface  easterlies  from  40  m/s  to 
30  m/s  occurs  due  primarily  to  the  dissipation  by  the 
momentum  flux  transfer.  The  maximum  easterly  jet  core  is 
raised  to  a  level  slightly  above  the  surface,  and  the 
strength  and  extent  of  the  narrow  lower  tropospheric 


westerly  region  appear  to  be  slightly  reduced  for  no-flux 
case . 

Near  the  surface,  the  cyclonic  vorticity 
decreases  in  the  with-flux  cases  (Figs.  3.12a  and  3.12c). 
This  indicates  that  the  weakening  of  the  easterly  maximum 
and  thermal  gradient  of  the  front  at  the  surface  does  lead 
to  a  corresponding  reduction  in  the  horizontal  wind  shear  in 
the  lower  troposphere. 

(2)  R . H . =  90%  Cases  (M90  vs.  M91) 

Similar  results  due  to  inclusion  of  the 
surface  fluxes  are  obtained  in  the  comparison  between  the 
with-flux  and  no-flux  cases  for  R.H.=90%  as  in  the  R.H.=60% 
cases : 

1)  The  difference  in  the  distribution  of  potential 
temperature  and  mixing  ratio  is  not  discernible, 
except  near  the  surface  to  the  north  of  front  (Figs. 
3.8b  and  3. 8d) ; 


2)  At  the 

surface, 

the 

inclusion 

of  moisture  and  heat 

fluxes 

weakens 

the 

meridional 

gradients 

of  mixing 

ratio 

(Figs. 

3.  Ob 

and  3 . 

8d)  and 

potential 

temperature  in 

the 

frontal  zone  (Figs. 

3.9b  and 

3. 9d) 

because 

of  the  stronger 

fluxes  in 

the  cooler 

and  drier  regions; 

3)  The  surface  momentum  transfer  leads  to  a  weakened 


surface  easterly  maximum  (Figs.  3.11b  and  3. lid), 
with  the  easterly  jet  core  occurring  slightly  above 
the  surf  ace ; 

4)  The  inclusion  of  surface  fluxes  causes  a  stronger 
and  more  concentrated  vertical  velocity  (Fig.  3.10) 
and  heating  rate  (Fig. 3. 14)  in  the  frontal  zone,  so 
that  the  upper  tropospheric  westerly  wind,  the 
westerly  regime  in  the  middle  and  lower  troposphere 
above  the  surface  front  are  stronger,  and  the 
meridional  wind  shear  in  the  lower  troposphere  is 
weaker  (Figs.  3.11b  and  3. lid  and  3.12b  and  3 . 1 2 d ) ; 
and 

5)  For  R.H.=903>,  the  multi-band  structure  in  vertical 
velocity  (Figs.  3.10b  and  3.10d),  for  the  with-flux 
case  is  less  conspicuous  than  that  for  the  no-flux 
case,  as  was  noted  earlier. 

C.  SUBTROPICAL  SIMULATIONS 

1 .  Typical  Subtropical  Cases  (S90) 

To  examine  the  time  evolution  of  subtropical 

experiments,  the  case  of  S90  (initial  R .  H  .  =  9  0  %  and  no 

surface  fluxes)  is  chosen  as  the  "typical  subtropical  case". 
The  initial  potential  temperature  and  mixing  ratio  are 
depicted  in  Fig.  3.16a  and  the  zonal  velocity  is  shown  in 

Fig.  3.16b.  The  choice  of  the  typical  subtropical  case  is 

motivated  by  an  intention  to  represent  the  major  differences 


between  the  mean  midlatitude  atmosphere  and  the  mean 
subtropical  atmosphere  during  the  Mei-Yu  season.  The  mean 
subtropical  Mei-Yu  environment  is  conditionally  unstable 
with  a  R.H.  around  85%.  Furthermore,  during  intense 
convection  periods  when  a  LLJ  is  usually  observed,  the  R.H. 
is  around  80%  or  higher  (Chen  and  Tsay,  1978).  While  these 
intense  episodes  are  local  and  transient,  and  therefore 
cannot  be  properly  studied  by  the  present  two-dimensional 
model,  they  are  the  only  periods  when  the  LLJ  is  clearly 
defined.  Therefore,  the  somewhat  higher  value  of  R.H. =90% 
is  chosen  so  that  its  results  may  also  shed  light  to  the 
situations  when  an  intense  LLJ  is  developed.  In  this  view 
the  "typical  subtropical  case"  may  also  be  considered  to 
represent  the  "typical  intense  LLJ  case." 

As  in  the  typical  midlatitude  case,  the  large-scale 
deformation  forcing  causes  an  accumulation  of  the  potential 
temperature  and  moisture  isolines  to  result  in  a  frontal 
structure.  However,  the  approach  towards  a  quasi-steady 
state  is  slower.  Furthermore,  all  fields  develop 
significant  small-scale  variations  on  the  equator  side  of 
the  front.  (For  example,  see'  the  potential  temperature  and 
mixing  ratio  fields  shown  in  Figs.  3.17a,b).  These 
variations,  which  are  similar  to  those  in  the  strong 
convection  cases  of  midlatitude  experiments,  are  apparently 
a  consequence  of  intense  cumulus  convection  in  the  model 
resulting  from  the  very  moist  and  unstable:  atmosphere. 


It  should  be  noted  that  once  the  initial  vertical 
0  e  profile  is  set  from  the  temperature  sounding  and 
relative  humidity,  it  will  remain  fixed  in  time  at  lateral 
boundaries.  Since  a  conditionally  unstable  profile  occurs 
in  the  lower  latitude  regions  (because  of  the  initial 
potential  temperature  and  mixing  ratio  distribution),  the 
instability  will  persist  on  the  equatorward  boundary 
regardless  of  convective  adjustment  in  the  interior.  This 
slower  development  into  a  quasi-steady  state  is  probably  due 
to  two  factors:  the  weaker  large-scale  deformation  forcing 
( D=0 . 5x10“ 5 s~ 1 )  prescribed  for  the  subtropical  cases  and  the 
presence  of  strong  conditional  instability  which  implies 
more  prominent  cumulus  convective  activity. 

The  small-scale  features  also  dominate  the  vertical 
velocity  field  such  that  the  frontal  structure  is  obscured 
by  high  frequency  oscillations.  To  remove  these  small-scale 
oscillations,  all  the  subtropical  experiments  fields  are 
subjected  to  a  daily  average  over  32  outputs  at  45  minute 
intervals.  Two  other  alternatives  were  considered:  to 
increase  the  control  of  large-scale  forcing,  or  to  decrease 
the  initial  relative  humidity.  However,  both  approaches 
will  reduce  the  contrasts  between  the  midlatitude  and  the 
subtropical  cases  especially  those  cases  with  clearly 
defined  LLJ,  which  is  one  of  the  main  objectives  of  this 
study.  Therefore,  it  was  decided  that  the  differences 
between  the  midlatitude  and  subtropical  conditions  should  be 


kept  as  large  as  possible  to  identify  the  underlying 
dynamics  of  the  development  of  the  LLJ.  Intermediate  cases 
with  R.H.=70%  and  80%  will  also  be  run  to  ascertain  the 
physical  significance  of  the  intensively  convective  R.H.=90% 
cases . 

The  daily  averages  of  the  potential  temperature  and 
mixing  ratio,  and  the  meridional  potential  temperature 
gradient  for  days  2-5  are  shown  in  Figs.  3.18  and  3.19, 
respectively.  Here  each  "day"  represents  the  average  of  the 
previous  24  hours  (i.e.,  day  5  is  the  average  between  96- 
120  h)  .  In  these  diagrams,  the  small-scale  features  have 
been  filtered  out,  and  the  development  of  the  frontal 
structure  is  clearly  evident.  The  surface  mixing  ratio 
gradient  and  range  agree  reasonably  well  with  Chen’s  (1978) 
observations  (Fig.  1.4b). 

The  meridional  potential  temperature  gradient  also 
evolves  into  a  frontal  structure  near  the  surface  around  16° 
after  4  days  (Fig.  3.19).  The  tilt  of  the  frontal  zone  is 
similar  to  the  inidlatitude  case  (Fig.  3.4),  but  the  well- 
defined  frontal  zone  is  confined  to  the  lower  troposphere. 
Chen  and  Chang’s  (1980)  comparison  between  different 
segments  of  an  elongated  Mei-Yu  front  that  extends  from  the 
subtropical  southern  China  to  the  midlatitude  central  Japan. 
They  observed  that  the  subtropical  segment  of  the  front  has 
less  vertical  tilt  and  is  shallower  than  the  midlatitude 
segment.  These  differences  were  attributed  to  the  stronger 


convection  in  the  subtropics.  It  will  be  shown  later  that 
in  the  subtropical  with-flux  case  (S91)  the  tilt  is  indeed 
less  while  the  vertical  extent  remains  shallow. 

By  day  2,  the  concentrated  rising  motion  (Fig.  3.20) 
develops  in  the  equatorial  latitudes  and  then  moves  poleward 
towards  the  area  of  maximum  deformation  forcing.  At  day  3, 
the  maximum  rising  motion  is  at  16°  with  a  magnitude 
~17xl0-4  cb/sec.  After  day  4,  it  shifts  to  18°  and  the 
maximum  decreases  to  4-5  x  10~4  cb/sec.  This  is  the  process 
for  the  model  to  develop  the  Mei-Yu  front  from  an  initial 
zero-perturbation  state.  In  the  real  atmosphere,  the 
process  is  much  slower,  where  the  monsoon  trough,  which  is 
developed  over  a  seasonal  time  scale,  is  already  in  place 
when  the  front  moves  into  the  Mei-Yu  position  from  the 
north.  The  development  of  a  slanted  structure  in  the  lower 
(z=0-3  km)  and  upper  parts  (5-10  km)  of  the  concentrated 
rising  region  is  evident.  The  slanting  is  northward  with 
height,  in  the  same  sense  as  the  slope  of  the  frontal 
surface.  The  multi-band  structure  is  not  clearly  seen  in 
the  24-h  averaged  vertical  velocity  field  but  is  more 
clearly  visible  when  instantaneous  (not  daily-averaged) 
vertical  velocity  fields  are  examined.  An  example  is  given 
in  Fig.  3.21,  which  is  the  instantaneous  vertical  velocity 
at  the  time-step  of  day  4.5. 

The  evolution  of  the  zonal  wind  is  illustrated  in 
Fig.  3.22.  Tile  gross  structure  of  an  upper  tropospheric 


westerly  jet,  poleward  of  a  lower  tropospheric  easterly  jet 
with  a  maximum  at  the  surface,  roughly  resembles  that  of  the 
typical  midlatitude  case  (M60,  Fig.  3.6).  However,  there 
are  significant  differences.  The  subtropical  westerly  jet 
core  is  lower  (at  10  km),  shallower  and  wider  in  the 
meridional  direction.  These  features  are  different  from 
what  would  be  expected  for  the  mean  polar  front  jet  and 
subtropical  jet.  Compared  to  Chen’s  (1978)  observation  of  a 
subtropical  Mei-Yu  case  (Fig.  1.4a),  the  jet  core  is  lower. 
The  reason  for  the  lower  jet  core  (at  10  km  versus  12  km  in 
Fig.  1.4a)  is  apparently  the  fixed  tropopause  which  is 
imposed  by  the  mean  sounding  at  10  km.  In  an  intense  Mei-Yu 
period,  the  deep  convection  and  the  tropopause  are  a  couple 
of  kilometers  higher,  where  the  jet  maximum  is  also  located. 
On  the  other  hand,  the  meridional  extent  of  the  simulated 
westerly  regime  and  that  observed  are  comparable.  A 
significant  easterly  regime  occupies  the  entire  tropical 
upper  troposphere  and  caps  a  westerly  regime  in  the  tropical 
lower  troposphere.  In  Fig.  1.4a,  the  observed  zonal  wind  is 
also  easterly  in  the  upper  level  tropics,  although  the  area 
is  somewhat  higher  and  more  equatorward  than  the  simulation. 
The  imposed  tropopause  in  the  model  is  again  the  reason  for 
this  difference. 

Within  the  tropical  low-level  westerly  regime  is  a 
remarkably  strong  and  narrow  jet  structure  centered  at  z=3 
km  at  the  latitude  of  maximum  rising  motion.  This  low  level 


jet  first  appears  at  day  2  with  a  magnitude  ~5  ra/s. 


It 


develops  to  >  10  m/s  at  day  5.  In  terms  of  the  relative 
position  with  respect  to  the  front,  this  low  level  jet  is 
similar  to  the  local  lower-tropospheric  westerly  regime 
simulated  in  the  midlatitude  cases,  but  the  intensity  is 
much  stronger.  This  intensity  agrees  reasonably  well  with 
the  observed  LLJ  strength  in  Fig.  1.4b,  especially  when  the 
addition  of  a  mean  westerly  is  considered. 

The  strong  horizontal  wind  shear  that  develops 
around  the  frontal  2one  is  shown  in  Fig.  3.23.  Compared  to 


Fig.  3.7,  this  subtropical  case  has  a  deep  lower 
tropospheric  shear  zone  that  extends  into  the  middle 

troposphere,  and  its  magnitude  poleward  of  the  low  level  jet 
significantly  exceeds  that  of  the  upper  tropospheric 

maximum.  The  slope  of  the  frontal  shear  zone  is  again 
i  larger  (less  vertical  tilt)  in  the  subtropical  case  than  in 

j  the  midlatitude  case  (Fig.  3.7),  at  least  in  the  lower 

troposphere.  In  Fig.  3.23d,  a  region  of  strong  negative 

relative  vorticity  exists  above  the  surface  front  near 

^  10  km,  which  indicates  the  possibility  of  inertial 

instability.  This  is  the  consequence  of  the  intense 

t  convection  and  may  not  be  physically  realistic.  It  is  not 

known  whether  this  instability  contributes  to  the  small- 
scale  variations. 

To  demonstrate  that  the  five-day  integration  gives  a 
reasonable  quasi-steady  state,  the  averages  of  each  field 
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over  the  first  five  days  of  integration  was  used  as  the 
initial  conditions  for  another  integration  of  120  h.  The 
results  of  potential  temperature  and  mixing  ratio, 
meridional  potential  temperature  gradient,  vertical  velocity 
and  zonal  wind  for  the  (daily  averaged)  new  day  5  are  shown 
in  Fig.  3.24a-d.  These  fields  are  similar  to  those  of  the 
original  daily  averaged  at  day  5  (Figs.  3.18,  3.19,  3.20  and 
3.22).  This  result  indicates  that  the  daily  averaged  fields 
are  stable  in  time,  in  spite  of  the  presence  of  small-scale 
variations  and  possible  inertial  instability. 

The  subtropical  frontal  positions  as  indicated  by 
the  various  fields  are  farther  equatorward  from  the  central 
axis  of  deformation  forcing  than  that  of  the  midlatitude 
cases.  In  the  midlatitude  case,  the  central  axis  of 
deformation  is  at  46°  and  the  frontal  structure  near  the 
surface  develops  at  about  42°.  In  the  subtropical  case,  the 
axis  of  deformation  is  at  22.5°,  while  the  surface  frontal 
zone  is  near  16°.  This  increased  distance  apparently 
reflects  the  effect  of  the  smaller  Coriolis  parameter  which 
leads  to  a  relatively  larger  ageostrophic  secondary 
circulation  around  the  frontal  zone.  The  advection  by  the 
lower  branch  of  this  enhanced  secondary  circulation  tends  to 
move  the  surface  front  farther  equatorward.  In  addition, 
the  effect  of  a  possible  secondary  circulation  south  of  the 


front,  that  is  induced  by  the  strong  convection,  may  also 


contribute  to  the  more  equatorward  position  of  the  surface 
f  ront . 

2 .  Other  Subtropical  Cases 

As  in  the  midlatitude  experiments,  the  effects  of 
varying  the  initial  relative  humidity  and  surface  fluxes  are 
studied  by  running  three  additional  subtropical  cases  (S60, 
S61  and  S91).  The  resultant  structure  at  day  5,  averaged 
over  the  preceding  24-h  period,  for  all  four  subtropical 
cases  are  shown  in  Figs.  3.25-3.29.  The  precipitation  and 
heating  rates  are  shown  in  Figs.  3.30-3.31.  In  these  cases, 
no  significant  small-scale  variations  are  present,  so  the 
averaging  is  unnecessary.  Nevertheless,  it  is  done  to  make 
the  analysis  procedure  of  all  subtropical  cases  consistent, 
a.  Effects  of  Moisture 

Panels  a  and  b  in  Figs.  3.25-3.29  are  for  S60 
(R.H.=60%,  no-flux)  and  S90  (R.H.=90S»,  no-flux)  cases, 
respectively.  The  R.H.=60%  case  (Fig.  3.25a)  has  a  less 
intense  moisture  front  throughout  the  lower  troposphere. 
The  poleward  vertical  tilt  of  the  maximum  mixing  ratio 

gradient  is  about  the  same  as  in  the  R.H.=90%  case  (Fig. 
3.25b).  However,  the  curvature  near  the  axis  of  the 

maximum  mixing  ratio  does  not  appear  as  sharp  as  in  the 
R.H.=903»  case  due  to  the  drier  air  and  weaker  gradient.  The 

short-  wave  feature  in  the  middle  troposphere  over  the 

surface  front  position  is  not  as  intense  as  in  the  R.H.=90% 
case,  but  it  is  still  very  notable,  which  indicates  the 
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presence  of  significant  cumulus  convection.  As  may  be 
expected,  the  R.H.=60%  case  has  a  lower  temperature  in  the 
upper  troposphere  (Figs.  3.25a  and  3.25b).  The  small-scale 
variations  near  the  surface  in  the  cold  air  (Fig.  3.25d- 
3.27d),  as  seen  in  midlatitude  cases,  evidence  the  shallow 
layer  of  conditional  instability  caused  by  surface  heat 
fluxes.  The  conditional  instability  can  be  examined  further 
by  calculating  the  0e  field  which  will  be  shown  later  in 
Fig.  4.6  (midlatitude  cases)  and  4.7  (subtropical  cases). 
The  weaker  meridional  potential  temperature  gradient  in  the 
frontal  zone  can  be  seen  by  comparing  Figs.  3.26a  and  3.26b, 
both  of  which  also  show  less  small-scale  features  elsewhere 
in  the  domain.  The  vertical  tilt  of  the  frontal  zone  is 
somewhat  larger  and  the  frontal  zone  extends  to  a  slightly 
lower  height  at  the  poleward  boundary. 

The  vertical  velocity  of  the  R.H.=60%  case  (Fig. 
3.27a)  has  a  concentrated  zone  of  rising  motion  in  the  lower 
and  middle  troposphere,  but  the  magnitude  is  less  than  in 
the  R.H.=90%  case. 

The  zonal  wind  of  the  R.H.=60%  case  is  much 
weaker  than  in  the  R.H.=90%  case  (Fig.  3.28a  versus  Fig. 
3.28b).  The  upper  tropospheric  westerly  jet  maximum  of 
<  25  m/s  is  only  about  70%  of  the  >  35  m/s  magnitude  of  the 
R.H.=90%  case.  This  is  consistent  with  the  expected 
difference  in  the  energy  conversion  between  the  two  cases, 
as  the  heating  rate  and  the  upward  motion  (which  correlates 


almost  completely  with  the  heatin)g  of  the  warm  air  is 
stronger  for  the  R.H.=902s  case.  The  former  is  also  situated 
more  equator ward.  The  low-level  westerly  jet  structure 
still  persists,  but  the  magnitude  is  reduced  by  half.  The 
upper  tropospheric  easterly  regime  almost  disappears,  and 
the  surface  easterly  maximum  is  only  reduced  slightly  (Fig. 
3.28a  and  Fig.  3.28b).  These  features  also  appear  in  the 
comparison  between  Figs.  3.29a  and  b,  which  primarily 
indicates  a  significant  reduction  of  upper  tropospheric 
horizontal  wind  shear. 

When  fluxes  are  included,  the  comparison  between 
the  R.H.=60%  (S61)  and  R.H.=90&  cases  (S91)  [panels  c  and  d 
in  Figs.  3.25-3.31,  respectively]  yields  similar  results. 
In  general,  the  meridional  temperature  and  mixing  ratio 
gradients  and  the  zonal  wind  are  stronger  in  the  R.H.=90% 
case . 

The  subtropical  R.H.=90%  cases  are  included  to 
study  the  intensely  convective  Mei-Yu  situations  when  the 
LLJ  is  fully  developed.  However,  the  high  initial  moisture 
and  the  subtropical  sounding  make  the  atmosphere  highly 
conditionally  unstable.  The  resultant  high  frequency,  small 
scale  variations  are  removed  by  the  24  h  averaging.  Even 
though  the  averages  remain  stable  for  extended  integrations, 
a  question  may  still  be  asked  as  to  whether  the  results  are 
not  distorted  seriously  by  the  small-scale  variations. 
Therefore,  two  intermediate  with-flux  cases,  R . H . =70%  (S71) 


and  R.H. =80%  (S81),  are  run  to  assess  the  physical 
significance  of  the  S91  case.  Fig.  3.33  compares  the  24-h 
averaged  zonal  wind  of  all  four  subtropical,  with-flux 
cases.  Here  it  can  be  seen  that  the  pattern  of  all  four 
cases  are  similar.  The  midlatitude  -upper  level  westerlies, 
lower  level  easterlies,  tropical  upper  easterlies  are 
present  in  all  cases,  and  the  intensity  of  these  features 
increases  systematically  as  the  initial  P.H.  increases  from 
60%  to  90%.  The  LLJ  is  indicated  by  a  very  small  area  of 
5  m/s  in  the  S61  case  (Fig.  3.33a),  this  area  expands 
considerably  with  increasing  humidity.  When  the  R.H. 
reaches  90%  (Fig.  3.33d),  the  jet  maximum  becomes  >  10  m/s. 
The  similarity  in  the  structure,  and  the  consistent 
intensification  of  the  circulation  with  increasing  moisture, 
indicate  that  the  subtropical  R.H. =90%  results  can  indeed  be 
viewed  as  physically  relevant  to  the  strong  LLJ  situations, 
b.  Effects  of  Surface  Fluxes 

When  surface  momentum,  heat  and  moisture  fluxes 
are  included  (Fig.  3.32),  these  important  contrasts  continue 
to  exist,  as  seen  by  comparing  panels  b  and  d  in  Figs.  3.25- 
3.31.  Similar  to  the  midlatitude  experiments,  the  addition 
of  surface  heat,  moisture  and  momentum  fluxes  leads  to  a 
higher  upper  tropospheric  temperature.  This  is  apparently 
due  to  the  increased  precipitation  and  the  associated  latent 
heat  release  (Figs.  3.30  and  3.31),  which  are  to  be  expected 


from  the  increased  sensible  and  latent  heat  supply  (Fig. 


3.32) .  A  weaker  potential  temperature  and  mixing  ratio 
gradient  at  the  surface  (Figs.  3.25  and  3.26)  is  due  to  the 
stronger  surface  fluxes  of  heat  and  moisture  in  the  cooler 
and  drier  region  (Figs.  3.32). 

The  vertical  velocity  of  the  with-flux  case 
displays  the  stacked  multi-band  structure  slanted  upward  and 
poleward  in  the  surface  frontal  zone  (Figs.  3.27b  and  d). 
Also,  the  magnitude  of  updraft  is  greater  than  in  the  no¬ 
flux  case.  This  is  consistent  with  the  precipitation  and 
heating  rates  (Fig's.  3.30  and  3.31)  and  is  due  to  the 
additional  moisture  supply  from  the  lower  boundary  (Fig. 

3.32)  . 

Comparisons  between  Figs.  3.28a  and  3.28c,  and 
between  Figs.  3.28b  and  3.28d  show  the  effects  of  surface 
fluxes  on  the  horizontal  circulation,  which  are  basically 
similar  to  the  results  of  the  raidlatitude  experiments.  When 
fluxes  are  included,  the  upper  tropospheric  zonal  winds, 
including  both  the  westerly  jet  and  the  less  intense 
easterly  jet,  are  enhanced.  The  lower  tropospheric  westerly 
jet  near  3-4  km  for  the  R.H.=90%  case  (Fig.  3.28d)  also 
intensifies.  Due  to  the  momentum  dissipation  at  the  surface 
(Fig.  3.32),  the  surface  easterlies  weaken,  and  the  jet  core 
is  shifted  slightly  upward  (Figs.  3.29b  and  d). 

The  comparisons  of  the  surface  fluxes  and  no 
fluxes  for  the  R  .  H  .  =  6  0  %  cases  (panels  a  and  c  of  Figs.  3.25- 
3.29)  basically  illustrate  similar  results. 


D.  SUMMARY  OF  THE  RESULTS 

The  results  of  the  foregoing  numerical  experiments  may 

be  summarized  as  follows: 

1)  In  the  midlatitude  case,  the  zonal  velocity  field  is 

dominated  by  a  westerly  regime  occupying  almost  the 
entire  upper  troposphere  and  an  easterly  regime 
occupying  the  lower  troposphere.  On  the  other  hand,  in 
the  subtropical  cases,  the  zonal  flow  equatorward  of  the 
front  tends  to  reverse,  with  expansive  easterlies  in  the 
upper  troposphere  and  weak  westerlies  in  the  lower 
troposphere.  A  small  regime  of  easterlies  develops  in 
the  midlatitude  upper  troposphere  only  in  the  high- 
humidity  cases. 

2)  At  the  latitude  of  the  maximum  rising  motion,  a  narrow 

westerly  wind  regime  occurs  in  the  lower  troposphere. 
In  midlatitude  cases,  it  appears  as  a  small  downward 

(and  equatorward)  penetration  of  the  upper  tropospheric 
westerlies.  In  the  subtropical  cases,  it  exists 

distinctly  as  a  well-defined,  westerly  maximum  similar 
to  the  low-level  jet. 

3)  The  poleward  tilts  of  all  fields  except  vertical 

velocities  are  more  prominent  in  the  midlatitude  cases 
are  either  comparable  or  stronger  than  in  the 
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subtropical  cases.  The  frontal  zone  also  extends  higher 
in  the  midlatitude  cases.  These  are  in  agreement  with 
observations  (e.g.,  Chen  and  Chang,  1980). 

The  inclusion  of  higher  initial  humidity  values  or  of 
surface-  fluxes  tends  to  give  stronger  vertical  velocity 
and  upper  tropospheric  zonal  winds.  This  is  a  result  of 
stronger  cumulus  convection  and  precipitation  whose 
latent  heat  release  provides  more  energy  input  to  the 
upper  troposphere.  At  lower  levels,  the  easteries  weaken 
as  a  consequence  of  the  momentum  dissipation  at  the 
surface . 

When  the  vertical  velocity  is  intense,  a  multi-band 
structure  often  appears.  It  is  characterized  by  more 
than  one  concentrated  updraft  region,  stacked  closely  to 
each  other  along  the  poleward  slanted  frontal  surface. 
The  stronger  updraft  is  situated  equatorward  and  has  a 
larger  vertical  extent. 

The  frontal  zone  develops  equatorward  of  the  dilation 
axis  of  the  large-scale  deformation  forcing.  In  the 
midlatitude  cases,  the  frontal  zone  is  closer  to  this 
axis  (about  4°  latitude  separation  at  the  surface), 
while  in  the  subtropical  cases  the  separation  is  greater 
(about  5-6°)  at  the  surface.  This  is  apparently  a  result 
of  the  smaller  f  and,  therefore,  a  relatively  larger 
ageostrophic  wind  which  advects  the  surface  front 


farther  equatorward. 


In  addition,  the  effects  of  the 


convection-induced  secondary  circulation  may  also 
contribute  to  the  more  equatorward  position  of  the 
surface  front. 

In  the  ensuing  section  these  numerical  results  will  be 
discussed  in  relationship  to  the  development  of  the  low- 


level  westerly  jet. 
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Fig.  3.1  Initial  potential  temperature  (solid),  8  (°K),  and  equivalent 
potential  temperature  (dash),  9e(°K),  soundings  at  the  center 
of  the  model  domain  for  midlatidude  cases  -  M60  and  M90 
and  subtropical  cases  -  S60  and  S90.  (see  text). 


Cross  section  of  mixing  ratio  (solid)  ,  q(g  kg‘1),  and  potential  temperature 
(dash),  0  (°K),  for  typical  midlatitude  case  M60,  at  (a)  day  1.0,  (b)  day  2.0,  (c) 


LATITUDE  •  LATITUDE 

Fig.  3.3,  except  for  horizontal  potential  temperature  gradient,  90/3 
-l).  Heavy  dash  lines  at  day  4,  (-0.02  °K  km"1)  indicate  frontal  zone. 
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Cross  section  of  mixing  ratio  (solid),  q(g  kg-i),  and  potential  temperature 
(dash),  8  (°K),  for  different  midlatitude  cases  at  day  4,  (a)  M60  MQr 


except  for  vertical  pressure  velocity,- (o(10'4  cb  s_1) 


except  for  vorticity,  5  (10'4  s_l). 
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Cross  section  of  24-hour  averaged  mixing  ratio  (solid),  q(g  kg_1)i  and  potential 
temperature  (dash),  0(°K),  for  (a)  day  2.0,  (b)  day  3.0,  (c)  day  4.0,  and  (d)  daj 


ere 


Fig.  3.15,  except  for  vertical  pressure  velocity,  -a)(10'4 


Cross  section  of  24-hour  averaged  mixing  ratio  (solid),  q (g  kr'1),  and  potential 
temperature  (dash),  0(°K),  for  different  subtropical  cases  (a)  S60,  (b)  S90,  (c) 
(d)  S91. 


except  for  horizontal  potential  temperature  gradient, 
i').  Heavy  dash  lines  (-0.01°K  km-1)  indicate  frontal  zone 


LATITU 


except  for  vorticity,  C(10-4 


Precipitation  rate  for  subtropical  cases  in  units  of  cm  day-*,  (a)  S60,  (b)  S90  (c)  S61  and 


Cross-section  of  total  heating  rate  for  subtropical  cases  in  units  of  degree 
S60,  (b)  S90  (c)  S61  and  <d)  S91. 


IV.  MECHANISM  FOR  THE  LOW  LEVEL  JET  ( LLJ ) 


The  previous  section  shows  that  a  low-level  westerly 
regime  tends  to  develop  in  the  frontal  zone  due  to  the 
intensifying  meridional  temperature  gradient  in  the 
lowest  few  kilometers.  In  many  of  the  subtropical  cases 
this  low-level  westerly  regime  actually  develops  into  a 
low-level  jet  whose  height  of  2-4  km  is  in  good  agreement 
with  the  levels  (800-600  mb)  of  the  observed  low-level 
jet  (Chen,  1977;  1978).  The  jet  maximum  of  10-15  m/s  in 
these  subtropical  cases  is  also  in  reasonable  agreement 
with  observations. 

In  the  Introduction,  it  was  pointed  out  that  there  are 
several  theories  for  the  development  of  the  low-level  jet 
(LLJ).  Since  effects  of  terrain  (Holton,  1967;  Wexler, 
1961)  and  diurnal  cycles  (Blackadar,  1957)  are  not 
included  in  the  model,  the  related  boundary-layer 


theories  cannot 

b  e 

treated  in 

this  study. 

Neither 

can 

the  theory  of 

the 

downward 

transport  of 

momentum 

by 

cumulus  convection  (Masoinoto  and  Ninomiya,  1969; 
Ninomiya,  1971),  because  the  convective  adjustment  used 
in  the  present  model  does  not  include  the 
parameterization  of  cumulus  momentum  transport. 

The  theory  proposed  by  Uccellini  and  Johnson  (1979)  in 
the  study  of  midlatitude  severe  storms  over  the  U.  S.  is 
based  on  a  wind-mass  adjustment,  whereby  the  secondary 
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circulation  forced  by  an  upper  level  jet  streak  causes  a 
vertical  coupling,  which  leads  to  the  development  of  a 
LLJ.  In  this  mechanism,  the  LLJ  will  develop  beneath  the 
upper  level  jet  streak  such  that  the  two  jets  will  be 
aligned  vertically  and  the  axes  will  be  perpendicular  to 
each  other.  Although  the  model-produced  low-level 
westerly  maximum  often  appears  (especially  in  the 
midlatitude  cases)  as  a  downward  extension  of  the  upper 
tropospheric  westerly  jet,  the  upper  level  jet  is  always 
situated  much  too  poleward  of  the  LLJ  core  for  wind-mass 
adjustment  mechanism  to  be  effective. 

In  the  present  study,  the  LLJ  in  the  subtropical 
experiments  is  much  better  defined  than  in  the 
midlatitude  experiments.  Furthermore,  in  all  cases  when 
the  vertical  motion  is  enhanced  due  to  increased  relative 
humidity  (therefore  also  enhanced  conditional 
instability)  and/or  surface  fluxes,  the  low-level 
westerlies  become  stronger.  This  is  also  true  with  the 
development  of  the  upper  level  easterlies  equatorward  of 
the  front,  which  also  appears  with  small  magnitudes  in 
the  midlatitude  experiments  with  higher  relative 
humidity.  In  all  subtropical  experiments,  these  upper- 
level  easterlies  expand  and  dominate  most  of  the  upper 
troposphere  eouatorward  of  the  front.  These  results 
suggest  that  when  the  vertical  motion  above  frontal  zone 
is  intense,  a  "reversed  Hadley  cell"  exists  equatorward 


of  the  front,  with  rising  motion  in  the  frontal  region 
and  sinking  motion  further  equatorward.  Such  an 
overturning  would  provide  an  easterly  Coriolis  torque  at 
the  upper-level  equatorward  flow  and  a  westerly  Coriolis 
torque  at  the  lower-level  poleward  flow,  which  may 
explain  the  upper-level  easterlies  and  the  low-level 
westerly  jet.  Indirect  evidence  of  the  "reversed  Hadley 
cell"  to  the  south  of  LLJ  may  be  found  in  an 
observational  study  by  Ninomiya  and  Akiyama  (1974).  They 
found  that  moist  air  is  located  above  the  LLJ  in  the 
middle  and  upper  troposphere,  but  the  air  south  of  the 
LLJ  is  fairly  dry  throughout  the  whole  troposphere.  This 
is  indicative  of  a  large-scale  subsidence  south  of  the 
front . 

More  evidence  of  this  "reversed  Hadley  cell"  is  seen 
in  the  temperature  distribution  of  all  cases  in  which 
strong  cumulus  convection  is  indicated.  For  example,  the 
temperature  in  Fig.  3.18  (S90  -  the  typical  subtropical 
case)  increases  equatorward  throughout  the  troposphere 
for  the  first  day.  Afterward,  reversed  meridional 
temperature  gradients  begin  to  appear  in  the  tropical 
region,  such  that  by  day  5  the  maximum  temperature  at  any 
level  in  the  middle-upper  troposphere  is  near  the  frontal 
zone.  This  situation  is  also  apparent  in  Fig.  3.25  for 
all  subtropical  cases.  The  reversed  temperature  gradient 
to  the  south  of  the  maximum  upward  motion  has  been 
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observed  both  in  the  Mei-Yu  front  (Chen  and  Chang,  1980) 
and  in  the  Baiu  front  (Matsumoto  et  ai.,  1970;  Ninomiya 
and  Akiyama,  1970).  Therefore,  the  equatorside  "reserved 
Hadley  cell"  is  actually  thermally  direct,  with  warm  air 
rising  in  the  frontal  zone  and  cool  air  sinking 
equatorward.  Hence,  this  circulation  can  be  maintained 
and  as  well  as  continue  to  provide  the  upper-level 

easterly  and  lowei — level  westerly  accelerations. 

The  next  question  concerns  the  identification  of  the 
main  mechanism(s)  for  this  equatorward  overturning. 

There  appears  to  be  three  possible  mechanisms: 

1)  The  secondary  circulation  normally  associated  with 
the  large-scale  flow  of  a  front; 

2)  The  cumulus  convection  due  to  the  conditional 

instability,  which  can  increase  the  latent  heat 

release  and  force  an  enhanced  secondary 

circulation; 

3)  The  third  possibility  is  based  on  the  observation 

that  the  vertical  velocity  field  of  several  of  the 
more  intense  cases  exhibits  a  multi-band 

structure.  Each’  of  the  rising  branches  is 

upright,  and  is  stacked  in  a  poleward  slanting 
configuration.  In  a  way  this  structure  resembles 
the  "slantwise  convection"  associated  with 
symmetric  instability  in  a  baroclinic  atmosphere 
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(e.g.,  Emanuel,  1985).  This  slantwise  convection 
may  also  be  a  relevant  mechanism. 

To  differentiate  the  last  two  mechanisms,  the  term 
"upright  convection"  will  be  used  to  describe  the 
basically  single-band,  vertically-oriented  convection 
which  is  normally  associated  with  conditional 
instability,  and  "slantwise  convection"  will  be  used  to 
describe  the  slanted  and  stacked  multi-band  structure 
which  resembles  the  convection  predicted  by  the  symmetric 
instability  theory.  Fig.  4.1  illustrates  schematically 
the  differences  among  the  three  mechanisms. 

The  meridional-vertical  streamfu notion  is  shown  in 
Figs.  4.2  and  4.3,  for  the  midlatitude  and  subtropical 
cases,  respectively.  For  the  midlatitude  R.H.=60%  cases 
(M60  and  M61),  the  atmosphere  remains  statically  stable, 
and  Figs.  4.2a,c  indicate  that  the  large-scale  frontal 
secondary  circulation  consists  of  a  single  cell  in  the 
Hadley  sense.  When  R .  H .  = 9  0  % ,  the  initial  sounding  of  the 
midlatitude  cases  (M90  and  M91  )  is  conditionally  unstable 
in  the  lower  latitudes.  At  day  4,  both  cases  show  a 
small  and  very  weak' "reversed  Hadley  cell"  in  the  middle 
troposphere  equatorward  of  the  front.  Apparently  this 
cell  is  responsible  for  the  small  region  of  upper  level 
easterlies  equatorward  of  the  front.  However,  the  modest 
cumulus  convection  through  the  release  of  conditional 
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instability  is  not  intense  enough  to  force  a  strong 
reversed  Hadley  cell  that  can  develop  a  strong  LLJ. 

All  the  subtropical  cases  are  initially  conditionally 
unstable.  A  larger  and  stronger  "reversed  Hadley  cell” 
is  present  in  all  cases  (Fig.  4.3).  This  cell  is 
situated  poleward  and  is  more  vertically  confined  for 
the  R  .  H  .  =  6  0  %  case  than  for  the  R  .  H  .  =  9  0  %  case  with  the 
middle- level  rising  motion  occurring  poleward  of  the 
surface  front.  This  is  the  situation  when  the  LLJ  is 
either  ill-defined  (S60  case,  Fig.  3.28a)  or  rather  weak 
(S61,  Fig.  28c).  When  R.H.=90&,  the  reversed  Hadley  cell 
is  located  more  equatorward,  is  more  intense  and  extends 
throughout  the  vertical  domain  (Figs.  4.3b,d).  The  LLJ 
is  well-defined  in  this  situation  (Figs.  3.28b,d),  and 
clearly  coincides  with  the  concentrated  poleward  flow  at 
the  lower  branch  of  the  reversed  Hadley  cell.  The 
cumulus  convection  is  the  strongest  in  these  cases  and  is 
apparently  responsible  for  generating  the  intense 
equatorside  meridional  circulation  which,  in  turn,  leads 
to  the  development  of  the  LLJ. 

It  is  clear  that  the  cross-frontal  secondary 
circulation  is  different  from  the  secondary  circulation 
equatorward  of  the  front,  and  therefore  it  can  not  induce 
the  LLJ.  However,  the  roles  of  upright  convection  and 
slantwise  convection  ere  not  easy  to  distinguish.  The 


slantwise  convection 


structure  is 


present  in  both 


midlatitude  (when  the  LLJ  is  not  present)  and  subtropical 
(when  the  LLJ  is  developed)  cases.  The  slanting  is 
apparently  related  to  the  structure  of  the  frontal 
surface.  Although  the  frontal  forcing  cannot,  by  itself, 
directly  cause  the  equatorward  secondary  circulation,  it 
exerts  a  strong  influence  on  the  environmental  conditions 
that  affect  convection.  This  may  involve  two  mechanisms. 
The  first  is  the  large-scale  lifting  of  moist  air.  This 
leads  to  either  stable  precipitation  or  a  "CISK" 
(Conditional  Instability  of  the  Second  Kind)  type 
convection,  where  the  low-level  convergence  is  provided 
by  the  frontal  forcing.  The  second  involves  the  creation 
of  conditions  for  symmetric  instability.  These  conditions 
occur  when  the  vertical  gradient  of  equivalent  potential 
temperature  (fie)  becomes  negative  along  the  surfaces  of 
pseudo-angular  momentum,  which  is  defined  by  the 
following  function  (Emanuel,  1983a, b): 

M  =  u  +  f  y . 

In  a  symmetrically  unstable  baroclinic  atmosphere, 
slantwise  convection  occurs  when  there  is  a  slantwise 
displacement  of  an  air  parcel  along  a  constant  M  surface. 
Based  on  a  diagnostic  semi-geos t rophic  model,  Emanuel 
(1985)  further  showed  that  even  in  the  presence  of  small 
symmetric  stability,  frontal  forcing  can  result  in  a 
strong,  concentrated,  sloping  updraft  which  is  located 


ahead  of  the  region  of  maximum  geostrophic  fronto- 
genetical  forcing.  This  possibility  is  supported  by 
Sanders  and  Bosart  (1985)  who  reported  observations  of  a 
concentrated  sloping  updraft  that  resembles  slant-wise 
convection  in  a  mesoscale  precipitation  band.  The  sloping 
updraft  remains  even  after  the  environment  becomes 
symmetrically  stablized.  Apparently  the  f ronto-genet ical 
process  can  force  the  moist  updraft  along  the  sloped  re¬ 
surfaces  where  narrow  convective  bands  can  develop. 
Thus,  slantwise  convection  may  be  expected  in  areas  of 
neutral  or  small  symmetric  stability  when  f rontogenet ical 
forcing  is  present. 

The  evolution  of  8e  and  the  M-function  for  the 
typical  midlatitude  and  subtropical  cases  (M60  and  S90) 
is  shown  in  Figs.  4.4  and  4.5.  In  the  midlatitude  case, 
the  f rontogene t ical  area  is  statically  stable  initially 
and  the  M-surfaces  tilt  only  slightly  (Fig.  4.4a).  As 
the  f rontogenet ical  process  proceeds,  the  frontal 
structure  develops  and  leads  to  a  stronger  tilt  of  the  M- 
surfaces  in  the  highly-baroclinic  frontal  zone.  This 
results  in  an  area  of  near-neutral  symmetric  instability 
after  day  1,  although  the  area  is  limited  to  the  lowest 
troposphere.  In  the  subtropical  case,  conditional 
instability  is  present  throughout  the  domain  initially 
(Fig.  4.5a).  This  condition  continues  until  day  4  when 
the  frontal  structure  is  developed  (Fig.  3.19)  and  the 


frontal  zone  becomes  stable  in  the  upright  convection 
sense.  As  a  result  of  f rontogenesis ,  the  M-surfaces  tilt 
more  towards  the  cold  air  in  the  frontal  zone.  This 
increased  tilt  leads  to  a  condition  of  symmetric 
instability  in  a  relatively  deep  layer  between  the 
surface  and  6  km  in  the  frontal  zone.  As  the  convective 
instability  is  eliminated,  a  slantwise  convection 
structure  develops  in  the  vertical  motion  field.  It 
appears  that  this  structure  is  produced  by  the 
frontogenesis  forcing  of  moist  updrafts  along  the  sloped 
M-surfaces,  similiar  to  Sanders  and  Bosart’s  (1985) 
observations . 

Distribution  of  day-5  8  <»  and  the  M-function  are 
presented  in  Figs.  4.6  and  4.7  for  the  midlatitude  and 
subtropical  experiments,  respectively.  In  the  R.H.=60% 
midlatitude  cases  (M60  and  M61,  Figs.  4.6a,c),  static 
stability  remains  positive  and  the  area  favorable  for 
slantwise  convection  is  rather  limited  and  is  confined  to 
the  lowest  2  km  in  the  frontal  zone.  Thus,  the  secondary 
circulation  is  almost  entirely  a  result  of  the  frontal 
forcing. 

Distribution  of  8e  and  the  M-function  are  given  in 
Figs.  4 . 6b , d  for  the  R.H.=90%  midlatitude  cases  (M90  and 
M91).  In  these  cases,  the  equatorward  side  of  the  front 
is  conditionally  unstable  while  the  poleward  side  is 


fro.ntal  zone  the  convective  adjustment  apparently  has 


stablized  the  atmosphere  in  the  vertical  direction. 
However,  the  9  *  and  M  lines  are  nearly  parallel  in  the 
frontal  zone,  and  the  area  favorable  for  sloping  updraft 
is  much  greater  than  that  for  the  R.H.=60&  cases.  In 
this  case,  a  low  level  maximum  wind  of  5  m/s  actually 
develops  near  z  =  3  km. 

Similar  to  the  R.H.=90%  midlatitude  cases,  all  the 
subtropical  without-flux  cases  (Fig.  4.7a,b)  have 
conditional  unstable  profiles  equatorward  of  the  front 
and  a  stable  0e  profile  poleward.  In  the  with-flux  cases 
(Fig.  4.7c,d),  the  conditional  instability  also  exists  in 
the  poleward  region,  a  result  of  the  supply  of  moisture 
and  heat  at  the  surface.  The  frontal  rising  motion  zones 

in  all  cases  again  have  stable  0  e  profiles,  but  shows 

signs  of  even  stronger  symmetric  instability  in  lower 
troposphere . 

In  the  S60  case  (Fig.  4.7a),  there  is  a  larger  area  of 
near  neutral  0e  profiles  along  the  M-surfaces  in  the 
vicinity  of  the  frontal  rising  motion.  In  the  S61  case 
(Fig.  4.7c),  the  9  e  profiles  in  the  lower  troposphere 
actually  become  unstable  along  the  M  surfaces  in  the 
rising  motion  zone.  The  same  is  true  for  the  S90  and  S91 

cases  (Figs.  4.7b,d),  except  the  area  of  symmetric 

instability  is  even  larger  than  in  the  R.H.=60%  cases. 
Thus,  for  all  cases  when  the  condition  of  symmetric 
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instability  is  suggested  by  the  6  e  profile  along  M- 
surfaces,  a  relatively  intense  LLJ  is  developed.  While 
in  the  present  model  framework,  it  is  difficult  to  assess 
directly  the  possible  roles  of  symmetric  instability,  the 
appearance  of  the  slantwise  convection  structure  seems  to 
be  an  important  factor  in  the  occurrence  of  the 


equatorward  secondary  circulation. 


35. 


24-hour  averaged  streamfunction,  >p  (mxcb  s-»)  at  day  5  for  different 


LATITUDE 


Cross  section  of  M-function  (solid)  in  unit  of  m  s_1  and 
equivalent  potential  temperature  (dash),  0e(°K)  for  typical 
midlatitude  case,  M60,  at  (a)  day  0.0,  (b)  day  1.0. 


(a) 


LATITUDE 

Fig.  4.5  Cross  section  of  24-hour  averaged  M-function  (solid)  in 
units  of  m  s'1  and  equivalent  potential  temperature  (dash), 
0e(°K)  for  typical  subtropical  case,  S90,  at  (a)  day  0., 
(b)  day  1.0. 
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55.0 


potential  temperature  (dash),  0  e(°K). 


SUMMARY  AND  CONCLUSIONS 


In  this  work  a  two-dimensional  frontal  model  is  used  to 
study  the  structure  and  behavior  of  the  Mei-Yu  front,  with 
special  emphasis  on  the  dynamics  of  the  low-level  jet  (LLJ) 
that  is  frequently  observed  ahead  of  the  front.  The  Mei-Yu 
front  develops  in  the  area  of  the  east  Asian  monsoon  trough 
during  early  summer,  and  is  the  most  important  rain- 
producing  system  over  east  Asia  outside  of  the  typhoon 
season.  This  study  is  motivated  by  the  unique  mixed 
midlatitude-tropical  properties  of  the  Mei-Yu  front  (e.g., 
Chen  and  Chang,  1980),  where  heavy  convection  and  strong 
low-level  horizontal  shear  are  associated  with  the  nearly 
zonal ly-ori ented  front.  The  occurrence  of  the  LLJ  is  highly 
correlated  with  heavy  convective  rainfall.  The  mechanism  for 
the  development  of  this  LLJ  is  the  main  focus  of  this  study. 

In  the  numerical  model,  the  time-varying  response  to  a 
constant  large-scale  forcing  was  treated  by  a  set  of 
perturbation  equations.  The  large-scale  forcing  was 
prescribed  by  a  stretching  deformation  field  which  is 
characteristic  of  the  f rontogenet ical  forcing  of  the  Mei-Yu 
front.  The  quasi-steady  state  responses  to  this  forcing 
were  obtained  by  integrating  the  perturbation  equations  from 
an  initial  state  which  is  given  by  the  seasonally-averaged 


zonal  flow.  Two  major  sets  of  experiments  were  conducted. 


T  h  e  first 


based  on  midlatitude  conditions,  where  the 


Coriolis  parameter  at  46°  latitude  and  an  initially  stable 
sounding  were  chosen.  The  second  is  based  on  subtropical 
conditions,  where  the  Coriolis  parameter  at  22.5°  latitude 
and  an  initial  sounding  that  is  conditionally  unstable  was 
chosen.  For  each  set  of  initial  conditions,  a  number  of 
supplemental  cases  were  integrated  to  study  the  effects  of 
varying  initial  humidity  values  and  surface  fluxes  of 
momentum,  heat  and  moisture.  In  the  midlatitude  cases,  the 
quasi-steady  state  is  reached  within  four  days  of 
integration.  In  the  subtropical  cases,  small-scale  features 
due  to  intense  convection  cause  the  motion  field  to 
fluctuate  throughout  the  integration.  A  24-hour  averaging 
was  applied  at  day  5  (96-120  h).  This  is  considered  the 
quasi-steady  state  as  further  integration  shows  that  the 
motion  remains  quasi-steady. 

For  both  midlatitude  and  subtropical  experiments,  the 
frontal  zone  develops  equatorward  of  the  dilation  axis.  The 
separation  between  the  front  and  the  dilation  axis  is 
greater  for  the  subtropical  cases  than  for  the  midlatitude 
cases.  This  larger  separation  is  apparently  due  to  the 
smaller-  f  and,  therefore,  a  relatively  large  ageostrophic 
secondary  circulation  around  the  frontal  zone  whose  lower 
branch  advects  the  surface  front  farther  equatorward.  The 
midlatitude  front  extends  deeply  into  the  upper  troposphere 
with  a  strong  poleward  tilt,  whereas  the  subtropical  front 
is  confined  to  the  lower  troposphere  with  less  tilt,  which 


is  in  good  agreement  with  observations  (e.g.,  Chen  and 
Chang,  1980).  Along  the  sloping  front,  slantwise  updrafts 
develop  with  a  multi-band  structure.  This  updraft  is  more 
evident  in  the  subtropical  cases  and  in  the  more  moist 
raidlatitude  cases. 

A  westerly  jet  in  the  upper  troposphere  and  an  easterly 
jet  in  the  lower  troposphere  develop  in  both  midlatitude  and 
subtropical  cases,  with  smaller  intensities  in  the 
subtropical  cases.  The  inclusion  of  higher  humidity  or 
surface  fluxes  leads  to  stronger  upper  westerlies, 
apparently  as  a  result  of  stronger  cumulus  convection.  The 
lower  level  easterlies  weaken  both  in  the  midlatitudes  and 
the  subtropics  with  the  inclusion  of  surface  fluxes.  This 
is  apparently  due  to  the  effect  of  momentum  dissipation  at 
the  surface.  An  interesting  feature  in  the  zonal  wind  field 
in  the  subtropical  cases  is  the  concurrent  development  of 
upper-level  easterlies  and  low-level  westerlies  equatorward 
of  the  front.  The  low-level  westerly  maximum  at  z  =  3-4  km 
r esemb  les  a  LLJ,  whose  intensity  increases  when  higher 
moisture  is  included.  In  the  midlatitude  cases,  only  a 
small  downward  and  equatorward  penetration  of  the  upper 
westerlies  appears  in  the  similar  area. 

While  several  theories  have  been  proposed  for  the 
development  of  low-level  jets  in  the  Mei-Yu  and  other 
circulation  systems,  the  choice  of  the  model  set  up  prevents 
consideration  of  some  of  these  theories  as  a  possible 


explanation  for  the  s imul at  ed  LLJ.  The  absence  of  terrain 


and  diurnal  radiation  cycles  in  the  model  precludes  the 
boundary- layer  theories  (Wexler,  1961;  Holton,  1967; 
Blackadar,  1957),  and  the  absence  of  a  cumulus  momentum 
transport  parameterization  precludes  the  mechanism  due  to 
cumulus  momentum  transport  (Masomoto  and  Ninomiya,  1969; 
Ninomiya,  1971).  In  addition,  the  vertical  coupling 
mechanism  suggested  by  Uccellini  and  Johnson  (1979)  also 
does  not  apply  in  the  Mei-Yu  simulation,  because  the  upper 
tropospheric  westerly  jet  is  located  too  far  from  the 
latitude  of  the  LLJ  (Fig.  1.4). 

The  appearance  of  an  easterly  regime  in  the  upper 
troposphere  on  the  equator  side  of  the  front  above  the  LLJ 
suggests  that  the  LLJ  may  be  the  result  of  a  secondary 
circulation  that  resembles  a  "reversed  Hadley"  cell.  This 
possibility  is  supported  by  the  temperature  distribution, 
which  indicates  that  the  "reversed  Hadley"  cell  is  actually 
a  thermally  direct  circulation  that  is  capable  of 
maintaining  the  kinetic  energy.  The  Lower  branch  of  this 
circulation  is  directed  poleward  toward  the  front,  and  it 
exerts  a  westerly  Coriolis  torque  which  can  develop  and 
maintain  the  LLJ.  A  meridional-vertical  s t rearaf unct i on 
reveals  the  structure  of  this  secondary  circulation,  with  a 
strong  lower  branch  return  flow  coinciding  with  the 
development  of  a  LLJ  in  the  more  moist,  subtropical  cases. 
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Several  possible  mechanisms  may  generate  the  reversed 
Hadley  cell.  The  possibility  that  the  reversed  Hadley  cell 
is  a  secondary  circulation  associated  with  the  basic  frontal 
structure  was  ruled  out  because  in  the  drier  midlatitude 
cases  the  frontal  secondary  circulation  is  a  single  cell 
around  the  front  with  upward  motion  on  the  equator  side  and 
downward  motion  on  the  pole  side.  No  LLJ  is  developed  in 
these  cases.  In  the  more  moist  (especially  the  subtropical) 
cases,  the  LLJ  developed  with  the  equator-side  secondary 
circulation.  This  suggests  the  importance  of  cumulus 
convection  in  developing  the  reversed  Hadley  cell  and  the 
LLJ. 

In  all  cases  where  a  LLJ  is  developed,  there  is  a 
sloping  updraft  with  multi-band  rising  motion  cells.  This 
suggests  that,  besides  regular  "upright"  convection  due  to 
conditional  instability,  slantwise  convection  arising  out  of 
symmetric  instability  may  also  be  relevant.  The  roles  of 
these  two  types  of  convection  were  studied  by  examining  the 
evolution  of  the  equivalent  potential  temperature 
distribution  in  the  vertical  and  along  the  pseudo-angular 
momentum  (M-function)  surfaces.  In  the  subtropical  cases 
where  the  atmosphere  is  conditionally  unstable  initially, 
upright  convection  develops  rapidly  and  the  conditional 
instability  is  quickly  released.  During  the  early  stages  of 


the  development,  the  M-function  surface  is  nearly  vertical, 
so  that  the  condition  for  symmetric  instability,  (38e/3z)M 


0,  is  also  given  by  that  for  conditional  instability,  (30e/ 
3  z )  y  <  0.  As  the  f rontogenet ical  process  continues,  the  M~ 
surfaces  in  the  frontal  zone  begin  to  slope  and  the 
condition  for  symmetric  instability  become  different  from 
that  for  conditional  instability.  As  the  conditionally 
unstable  soundings  adjusts  to  neutral,  the  0e  sounding 
along  M  surfaces  is  either  neutral  or  slightly  unstable  in 
the  middle-lower  troposphere  below  6  km.  At  this  stage, 
slantwise  convection  is  observed  in  the  frontal  zone  and  is 
indicative  of  the  possible  importance  of  symmetric 
instability  processes. 

In  the  more  moist  midlatitude  cases,  slantwise 
convection  is  also  observed.  In  these  cases,  the  frontal 
zone  remains  statically  stable.  However,  the  f rontogenet ic 
forcing  produces  slanted  M  surfaces  there,  with  the  0? 
profiles  on  these  surfaces  reduced  to  only  slightly  stable 
based  on  the  symmetric  instability  criterion.  This  result  is 
in  agreement  with  the  recent  theory  of  Emanuel  (1985)  who 
showed  that,  under  f r on t ogene t i c  forcing,  sloping  updrafts 
can  occur  in  a  concentrated  region  on  the  warm  side  of  the 
frontal  zone,  and  result  in  an  enhanced  secondary 
circulation.  It  also  agrees  with  the  observation  by  Sanders 
and  Bosart  (1985)  of  a  mesoscale  precipitation  band  in  a 
winter-time  extratropical  cyclone  over  northeastern  United 


States,  in  which  a  sloping  updraft  was  observed  in  an 


environment  in  which  the  symmetric  stability  was  weakly 
pos i t i ve . 

The  main  findings  of  this  study  are  that  the  LLJ  in  a 
subtropical  Mei-Yu  front  can  develop  as  a  result  of  the 
Coriolis  torque  of  a  secondary  circulation  on  the  equator 
side  of  the  front.  This  secondary  circulation  is  thermally 
direct  and  is  in  the  reversed  Hadley  sense.  It  is  driven  by 
intense  condensation  heating  due  to  heavy  precipitation. 
The  upward  motion  takes  the  form  of  slantwise  convection  as 
a  result  of  the  fronto genesis  forcing.  The  slantwise 

convection  begins  to  appear  after  the  atmosphere  is 
stablized  by  the  convective  adjustment.  It  is  conceivable 
that,  without  the  slantwise  convection,  which  is 
concentrated  along  the  frontal  surface,  the  latent  heat 
release  would  spread  out  over  a  broader  area.  This  would 

give  rise  to  a  weaker  secondary  circulation.  These 
conclusions  are  consistent  with  the  observations  that  low- 
level  jets  are  often  associated  with  intense  convection  and 
heavy  rainfall  in  the  Mei-Yu  front,  and  the  indirect 
evidence  of  a  sinking  region  south  of  the  Baiu  front  as 
revealed  by  Mats  urn oto’s  (1972)  moisture  analysis.  Further 
study  is  needed  to  clarify  the  possible  roles  of  the 
symmetric  instability  mechanism.  However,  the  strong  in- 

phase  correlation  between  the  rising  velocity  and  heating 
indicates  that  the  main  conclusion  concerning  the  role  of 
the  secondary  circulation  should  not  depend  on  the  details 


of  the  formulation  of  condensation  heating  in  the  model.  In 
fact,  the  frontal  forcing  is  the  dominating  factor  that 
determines  the  region  of  maximum  rising  motion.  Given  the 
moist  subtropical  environment  for  the  Mei-Yu  system,  if  a 
specified  (instead  of  predicted)  heating  function  is 
appropriately  prescribed  in  the  model,  a  strong  secondary 
circulation  equatorward  of  the  front  may  also  be  expected. 

Observational  studies  (e.g.,  Chen  and  Tsay,  1978) 
reported  that  the  Mei-Yu  convection  and  heavy  rainfall  are 
organized  in  mesoscale  convective  complexes  (MCC),  which  can 
not  be  treated  in  the  present  two-dimensional  model.  A 
three-dimensional  mesoscale  numerical  predication  model, 
where  the  MCC’s  can  be  explicitly  simulated,  should  be  used 
to  reexamine  the  relationship  between  the  LLJ,  secondary 
circulation  and  convection.  A  three-dimensional  model  is 
also  necessary  to  study  the  entire  life  cycle  of  the  Mei-Yu 
front  and  the  LLJ.  The  development  and  decaying  stages  of 
the  LLJ  and  the  Mei-Yu  front  are  of  different  time  scales. 
While  the  Mei-Yu  front  usually  persists  for  several  days  to 
one  week  or  more,  the  appearance  of  the  LLJ  is  more 
transient  in  nature.  The  time  evolution  of  the  LLJ  within 
the  quasi-steady  Mei-Yu  front  is  therefore  an  interesting 
topic  for  further  studies. 

The  model  may  also  be  improved  to  study  additional 
effects.  The  inclusion  of  a  planetary  boundary  layer 
parameterization  that  is  more  sophisticated  than  the  bulk 


aerodynamics  and  is  coupled'  with  a  diurnal  radiation  cycle 
will  better  represent  the  effects  of  friction,  especially 
the  interaction  with  radiation  processes.  The  inclusion  of 
terrain  may  be  particularly  important,-  since  the  Mei-Yu 
front  is  over  the  mountainous  region  of  southern  China, 
Taiwan  and  Japan.  The  topographical  effect  of  the  mountains 
and  the  land-sea  complex  in  this  region  undoubtedly  play  an 
important  role  in  both  the  large-scale  front  and  the 
embedded  mesoscale  convective  systems.  In  addition,  Keyser 
and  Anthes  (1982,  1986)  have  argued  that  a  stratification- 
dependent  boundary  layer  effect  is  important  for  the 
intensity  of  vertical  velocity  in  a  surface  front,  while  Mak 
and  Bannon  (1984)  argued  that  condensation  heating  alone 
plays  a  more  critical  role.  A  sophisticated  planetary 
boundary  layer  parameterization  may  help  to  clarify  this 


issue.  Finally,  further  observational  studies  are  also 
necessary  to  investigate  the  Mei-Yu  processes  both  in  the 


Appendix 


A.  SPECIFICATION  OF  MODEL  GRIDS 

The  horizontal  grids  are  uniform.  The  grid  interval  is 
60  km.  The  pressure,  p,  is  used  as  vertical  coordinates  and 
k  is  the  vertical  grid  index.  The  pressure  of  A"th  level, 
Pk  ,  in  the  model  is  specified  by  the  exponential 
relationship  with  height,  Zk , 

Pk  -  P„  e'Zk  (A.l) 


where  po  is  the  pressure  at  the  model  bottom. 

Apply  (A.l)  to  the  total  height  of  the  model,  zt , 


where  pt  is  the  pressure  at  model  top. 

Using  uniform  increment  of  height  in  vertical, 
corresponding  pressure  is 


P 


k 


(k-1) 


Jin 


the 


where  kt  is  the  total  number  of  vertical  levels  in  the 


model . 


In  this  study,  the  total  number  of  vertical  levels,  kt  , 


is  54 . 

The 

corresponding  pressure  levels  used 

i  n 

the  mode  1 

are  shown  in 

Table  A. 

1 

TABLE 

A .  1  FORECAST' 

LEVELS  ( cb ) 

LEVEL 

p(cb) 

P  (Cb) 

P  (cb) 

1 

95.86 

19 

44.81 

37 

20.95 

2 

91.9 

20 

4  2.96 

38 

20.08 

3 

88.1 

21 

41.18 

39 

19.25 

4 

84.45 

22 

39.48 

40 

18.25 

5 

80.96 

23 

37.85 

41 

17.69 

6 

77.61 

24 

36.28 

42 

16.96 

7 

74.4 

25 

34 . 78 

43 

16.26 

8 

71.32 

26 

33.34 

44 

15.59 

9 

68.37 

27 

31.96 

45 

14.94 

10 

65.54 

28 

30.64 

46 

14.32 

11 

62.83 

29 

29.37 

47 

13.73 

12 

60.23 

30 

28  . 16 

48 

13.16 

13 

57.54 

31 

26.99 

49 

12.62 

14 

55.35 

32 

25.88 

50 

12.10 

15 

53.06 

33 

24.81 

51 

11 . 60 

16 

50.87 

34 

23.78 

52 

11.12 

17 

48.76 

35 

22.80 

53 

10.66 

18 

46.75 

36 

21.85 

54 

10.22 
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C.  FINITE  DIFFERENCE  FORMULATION 

The  equation  set  (2.1.16)  and  ( 2 . 1 . 19 ) - ( 2 . 1 . 23 )  can  be 
integrated  by  a  pure  inarching  process.  In  this  study  a 
spatially  staggered  finite-difference  scheme  is  used  as 
compared  to  the  unstaggered  scheme  used  in  Williams  (1974). 
The  variables  6,  it  and  q  are  carried  at  the  same  points  in 
y,  and  u  and  v  are  carried  equidistant  between  these  points 
(see  Fig.  B.l).  Winninghoff  (1968),  Arakawa  and  Lamb  (1977) 
and  Schoenstadt  (1980)  have  shown  that  this  staggered  grid, 
which  is  known  as  scheme  B,  is  much  better  than  the 
unstaggered  grid  (scheme  A)  for  geostrophic  adjustment. 
This  scheme  is  especially  superior  when  small-scale  heating 
occurs,  such  as  near  a  front.  The  nonlinear  terms  are 
arranged  to  conserve  mean  squares  in  the  advection  terms. 
Leapfrog  time  differencing  is  used,  except  the  model  is 
restarted  every  fifth  timestep  with  an  Euler  backward  step 
following  Arakawa  and  Lamb  (1977). 

For  the  finite  difference  formulation,  we  discuss  only 
the  spatial  derivative  terms  in  the  u  and  v  equations  [i.e., 
(2. 1.16)  and  (2.1.19)], 

^T7vF]j,'<  =  [(vJ+l,k  +  vJ,k)  (FJ+1  +  Fj,k) 


•  (vJ,k  +  vj-l »k>(Fj ,k  +  Fj-l,k^/4(yj+l  -  yj) 


CSy“Flj,k  *  fCj+l.Ktl  +»j>k+1)  (Fjlk+1  +  Fjjk) 

-  +  “j,k>  (fj,k  ♦  Fj,k-l)V4(pk+1  -  Pk) 


=  f F j *1 , k  -  Fj-I,k)/Gj+1  *  yj-i) 


■  t  -  (Fj-k  - 


fyj+1  -  yj) 


-  yj-1) 


^4/(yj+l  -  yj) 


[4-  K  ^1,  k  .  r  <M  (FJ.^  -  ri.O 

ap  ap-'J.*  ~  ~r~ 

(Pk+l  -  Pk) 


K^Fj,k  "  Fj,k-l) 
(Pk  ~  Pk-1 ) 


( Pk+l 


where  F  is  either  u  or  v 


Pk) 
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For  the  9  and  q  equations, 


vG:1j,l<  B  Cvj,k(Gj+l,k  +  Gj,k) 


vj-i »k^Gj ,k  +  Gj-i,k)/2(yj  -  yj-i) 


«G]j 


,k 


C(“jtk+l(fij,lc+l  +  Gj ,0  ”^j,k(Gj»k  +  Gj,k-l)/2(Pk+l  "  Pk) 


^j.k  =  fRj+i,k  -  Gj-i,k)/(yj+i  +  yj-i> 


(Gj+l,k  -  Gj,k) 

‘  (yj+i-yj) 


(Gj,k  “  Gj>l,k) 

-(y.  )~ 


•]/(yj  -  yj-i) 


dp 


Kffti.k  *  t- 


Kk+1  ^Gj,k+1  "  Gj,k)  Kk(Gj,k  *  Gj,k-l) 


(Pk+1  "  Pk) 


(Pk  '  Pk -1 ) 


]/(Pk+l  '  Pk) 


where  G  is  either  0  or  q. 
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